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HOW BAD/GOOD ARE THE EXTERNAL FORWARD SHOCK 
AFTERGLOW MODELS OF GAMMA-RAY BURSTS? 
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Deng 1,3 , Song-Mei Qin 7,2 , Qing-Wen Tang 8 , D. Alexander Kami 9 , Felix Ryde 10 , Pawan 

Kumar 11 

ABSTRACT 

The external forward shock models have been the standard paradigm to inter¬ 
pret the broad-band afterglow data of gamma-ray bursts (GRBs). One prediction 
of the models is that some afterglow temporal breaks at different energy bands 
should be achromatic , namely, the break times should be the same in different 
frequencies. Multi-wavelength observations in the Swift era have revealed chro¬ 
matic afterglow behaviors at least in some GRBs, casting doubts on the external 
forward shock origin of GRB afterglows. In this paper, using a large sample of 
GRBs with both X-ray and optical afterglow data, we perform a systematic study 
to address the question: how bad/good are the external forward shock models? 
Our sample includes 85 GRBs up to March 2014 with well-monitored X-ray and 
optical lightcurves. Based on how well the data abide by the external forward 
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shock models, we categorize them into five grades and three samples. The first 
two grades (Grade I and II) include 45/85 GRBs. They show evidence of, or 
are consistent with having, an achromatic break. The temporal/spectral behav¬ 
iors in each afterglow segment are consistent with the predictions (the “closure 
relations”) of the forward shock models. These GRBs are included in the Gold 
sample. The next two grades (Grade III and IV) include 37/85 GRBs. They 
are also consistent with having an achromatic break, even though one or more 
afterglow segments do not comply with the closure relations. These GRBs are 
included in the Silver sample. Finally, Grade V (3/85) shows direct evidence of 
chromatic behaviors, suggesting that the external shock models are inconsistent 
with the data. These are included in the Bad sample. We further perform sta¬ 
tistical analyses of various observational properties (temporal index a, spectral 
index /3, break time R) and model parameters (energy injection index q , electron 
spectral index p, jet opening angle 6j, radiative efficiency r^ 7 , etc) of the GRBs 
in the Gold Sample, and derive constraints on the magnetization parameter cb 
in the forward shock. Overall, we conclude that the simplest external forward 
shock models can account for the multi-wavelength afterglow data of at least half 
of the GRBs. When more advanced modeling (e.g., long-lasting reverse shock, 
structured jets, arbitrary circumburst medium density profile) is invoked, up to 
> 90% of the afterglows may be interpreted within the framework of the external 
shock models. 

Subject headings: radiation mechanisms: non-thermal — gamma-rays: bursts - 
method: statistics 


Introduction 


Gamma-ray bursts (GRBs) are the most luminous explosions in the universe. They 
signify the birth of a stellar-mass black hole or a rapidly rotating magnetized neutron star 


during core collapses of massive stars or mergers of compact objects (Kumar & Zhang 2015 
for a recent review). 


Multi-wave length GRB afterglows were predicted (IMeszaros fe Rees 


first discoveries f Costa et al. 1997 : van Paradiis et al. 19971 : 


Frail et ah 


99.71) before their 


1997 ). This was based 


on a generic external forward shock model. Regardless of the physical nature of progenitor 
and central engine, a relativistic jet is launched, which is decelerated by a circumburst 
medium by a pair of external (forward and reverse) shocks. The reverse shock is likely short¬ 
lived. The forward shock, on the other hand, continues to plough into the medium as the 
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jet is decelerated. Synchrotron radiation of electrons accelerated from the external forward 
shock powers broad-band electromagnetic radiation with a decreasing amplitude. This is 


the broad-band afterglow of GRBs ( 

Meszaros & Ree 

s 1997; 

Sari et al. 

1998; 

Meszaros et al. 

1998; 

Rhoads 

1999; 

Sari et al. 

1999; 

Chevalier & Li 

20001). 


Before 2004, the observations of the broad band late time afterglow emission of GRBs 


ies (e.g., 

Wiiers et al. 

1997; 

Waxrnan 

1997 

. 

vViiers & Ga 

ama 

1999; 

Harrison et al. 

1999 

Huang et al. 

4999, 

2000J; 

Panaitescu & Kumar 

2001, 

2002; 

Yost et al. 

2003; 

Wu et al. 

2004) 


suggested that these late-time data are generally consistent with the predictions of the ex¬ 
ternal forward shock models. 


The launch of the Swift, satellite in 2004 (Gehrels et ah| 2004) allowed systematic observa¬ 
tions of the multi-wavelength GRB afterglow at early epochs. These data, especially the early 
X-ray afterglow data, presented surp rises to modelers. The overall X-ray lightcurves include 


five distinct temporal components ( Zhang et al. 


conne cted to the prompt emission ( Tagliaferri et al. 2005 : Barthelmv et al. 112005 : Zhang et al. 


2006jh I: an early time steep decay phase 


2007d l: II: a shallow decay (or pl ateau) phase, which may signify continuous energy injec - 
tion of energy into the blastwave fjZhang et al.ll2006l; iNousek et al.ll2006t iLiang et al.l 120071) : 
III: a normal decay phase consistent with the forward shock emission of a constant-e nergy 


^ ^ — 1 

Racusin et al. 

2009): and V: erratic X-ray f 

ares, likely powered by late ce 

ntral engine activ- 

ities 

floka et a 

1. 20^ 

Burrows et al. 2005c 

Fan & Wei 2005]; Zhang et al. 

2006s; Liang et al. 

2006; 

Lazzati & Perna 

2007: Chincarini et al. 2007: Maxham & Zhang 2009: Margutti et al. 

2010) 

. The components I and V are are believed to be of an internal origin (in contrast to 


the external shock origin). The other three components (II, III and IV) may be interpreted 
within the framework of the external shock models. 


The optical afterglow light curves also show in 

terest 

ing temporal behaviors ( 

Liang et al. 

2006; Nardini et al. 

2006; 

kauri et al. 2006. 2010, 

2011 

Panaitescu & Vestrand 

20oi 2011; 

Li et al. 2012; 

Liang et al. 

2013 

Wang et al. 2013 

; Yi et al. 2013). In the simi 

iar spirit as 
ght curveof 

Zhang et al. ( 

2006). 

Li et al. (2 

012) attempted to summarize a “synthetic” li 


optical emission. They found more components with distinct physical origins: la: prompt 
optical flares; lb: an early optical flare of an external reverse shock origin; II: an early 
shallow-decay segment; III: the standard afterglow component (the normal decay component, 
sometimes with an early onset hump); IV: the post-jet-break phase; V: late optical flares; 
VI: late rebrightening humps; and VII: late supernova (SN) bumps. The components II, III 
and IV can find their counterparts in the canonical X-ray light curve (components II, III, 
and IV in IZhang et ah 2006J). Some flares in the optical band have counterparts in X-rays, 
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but some others do not (ISwenson et ah 20131 ). Some components (e.g., the reverse shock 
component lb and the supernova component VII) are unique for the optical band only. 

There are two types of temporal breaks in the external shock models. One t ype co r¬ 


respo nds to the crossing of a characteristc frequency in the observational band flSari et ah 


19981 1. Such spectrally-related breaks occur at different epochs in different energy bands, and 
therefore are chromatic. A testable feature of such a break is that the spectral indices before 
and after the temporal break should be distinctly different. The second type of breaks are 
related to the hydrodynamic or geometric properties of the system. Since both effects affect 
the global behavior of the blastwave, these breaks should be achromatic, i.e. the temporal 
breaks in different energy bands should occur around the same observational time. 

Most observed breaks in the GRB lightcurves are likely of a hydrodynamic or geometric 
origin. Observationally, essentially all the temporal breaks observed in the X-ray lightcurves 
are consistent with having no spectral changes across the break times flLiang et al.l 2007, 
20081 ). Theoretically, the spectral breaks, especially the cooling break, are predicted to 
be very smooth, and are barely observable fro m the data ( IThm V Zhang 2014al . see also 
Granot fe Sarill2002l; Ivan Eerten fe Wii era 120091) . As a result, one expects that the temporal 
breaks seen by Swift should be strictly achromatic based on the external forward shock 
models. 

Broad-band afterglow data of GRBs are rapidly accumulating. Shortly after Swift 
detected early X-ray afterglow lightcurves of GRBs, some authors noticed that the ba¬ 
sic r equirement of achromaticity of GRB afterglows is violated at least in some GRB s 
(e.g., Panaitescu et al. 2006 : Fan fe Piran 2006 : Huang et al. 2007 : Liang et al. 2007 . 2008 1. 
In particular, while a significant break is seen in the X-ray lightcurves of some GRBs, 
the optical lightcurve does not show evidence of a break at the corresponding time (e.g., 
Troia et al.l 2007 1: Molinari et ah 20071 ). Such a puzzling effect led theorists to suggest vari¬ 
ous non-forward-shock models of the X-ray a fterglow: the long-lasting reverse shock model 
( Genet et al. 2007 : Uhm &; Beloborodov 2007), the dust scattering model (Shao fc Dai 2007), 
and the long-lasting central engine model ( Ghisellini et al. 2007 : Kumar et al. 2008a b|). 


Indeed, if most GRB afterglows are chromatic, one must throw away the standard 
forward shock paradigm, and probably attribute other factors, in partular, the long-lasting 
central engine, to account for the X-ray afterglow. This would have profound implications 
for our understanding of the GRB central engine and emission physics. Yet, there seem 
to exist some GRBs (e.g., the latest bright GRB 130427A) whose multi- wave length data 
are consistent wit h the simplest forward shock afterglow model (e.g., iMaselli et al.l I20l4 
Perlcv et al. 2014 ). 
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It is therefore natural to ask the following question: in general how bad or how good 
are the external forward shock models in interpreting the GRB afterglow data? 


This paper aims at addressing this question through a systematic data analysis and 
theoretical modeling of a large sample of multi-wavelength afterglows. We study a sample 
of 85 Swift GRBs up to March 2014, which all have high-quality X-ray and optical light 
curve data to allow us to study the compliance of the data to the external forward shock 
models. The sample selection and data analyses are described in §2. The theoretical external 
forward shock model, in particular, the so-called closure relations are presented in §3. In 
§4, we grade the afterglows based on how well they abide by the forward shock models, 
and categorize them into five grades and three samples. A statistical analysis of various 
observational and theoretical parameters for the Gold sample is presented in §5. Our results 
are summarized in §6 with some discussion. We notice that Li et al. 02015 ) recently carried 
out a similar analysis, with the focus on the consistency of the data with afterglow models in 
individual temporal segments of X-ray and optical lightcurves, without analyzing the global 
achromatic/chromatic behaviors of the afterglows. 


Throughout the paper, the subscripts “O” and “X” denote the optical and X-ray band, 
respectively, and the subsripts “1” and “2” denote the pre- and post-break segments, respec¬ 
tively. In addition, two spectral regimes are defined: “I” for v > max(y m , v c ), and “IF' for 
I'm < v < u c , where v m and v c are the minimum injection frequency and cooling frequency 
for synchrotron radiation, respectively. 


2. Sample and Data 

We systematically investigate all the Swift GRBs that have X-ray and optical afterglow 
data, over a span of almost 10 years from the launch of Swift to March 2014. A sample 
of ~260 optical light curves are compiled from published papers or GCN Circulars, and a 
sample of ~900 X-ray light curves are obtained from the Swift XRT data archive. Well- 
sampled light curves in both X-ray and optical bands are available for 99 GRBs. Fifteen 
GRBs do not have well constrained spectral indices either in optical or in X-ray bands to allow 
us to perform some theoretical constraints (see details below). Fourteen of them are removed 
from the sample. GRB 070420 is the only GRB without adequate spectral information that 
is included in our sample. This is because it has a clear chromatic feature, which allows us 
to group it into the Bad sample even if the spectral information is not available (see details 
in 44.2ft . The remaining 84 GRBs are included in our final sample, whose information is 
presented in Table [0 For the optical data, the correction due to Galactic extinction is taken 
into account using the reddening map presented by Schlegel et al. (1998). Due to large 
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uncertainties, we do not make corrections to the extinction in the GRB host galaxies. 


In order to quantify the rich temporal features of GRB lightcurves, we fit the lightcurves 
with a model of multiple components. The basic component of our model is either a single 
power-law (SPL) function 

Fx = F m t~ a (1) 

or a smooth broken power-law (BPL) function 


F'2 — F 02 



( 2 ) 


where a, an, a 2 are the temporal slopes, R is the break time, and u measures the sharpness 
of the break. In some afterglow models, a double broken power-law light curve is expected. 
For example, it is theoretically expected that the afterglow light curve may have a shallow 
segment early on due to energy injection, then transits to a normal decay segment when 
energy injection is over, and finally steepens due to a jet break (e.g., in the canonical X-ray 
afterglow lightcurve, Zhang et al. 2006). We therefore also consider a smooth triple-power- 
law (TPL) function to fit some lightcurves. In these cases, we extend equation (|2]) (with t \, 
defined as Gy) to the following function flLiang et ah| 2008) 


F 3 = (F^ 2 + F^ 2 )~ 1/U2 


(3) 


where cu 2 is the sharpness factor of the second break at G, 2 5 and 

Ft = F 2 (t h _ 2 ) (J-) 3 . (4) 

We perform best fits to the data using a subroutine called MPFIlQ. The sharpness 
parameter u is usually adopted as 3 or 1 in our fitting. The parameter G is not significantly 



One focus of our analysis is to study the “chromaticity” of the lightcurves in the X-ray 
and optical bands. In principle there are two approaches to do this. The first approach is 


1 http://www.physics.wisc.edu/ craigm/idl/fitting.html. 
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to blindly search for tb using the best fits to the optical and X-ray data, respectively, and 
compare how different t 
times in the two bands 


Liang et al. 

2007, 

2008; 

Li et al. 

2012. 

2015) 


is to start with the achromatic assumption and investigate how bad the data violate such 
an assumption. By doing so, we reduce one free parameter, and impose a same tb in both 
bands in the model. We believe that this second approach is more reasonable to address the 
question “how bad the external forward shock models are”, so we adopt the second approach 
with the assist of the first approach. The detailed procedure of our light curve fitting is as 
follows: 


• For each GRB, we first fit the optical and X-ray afterglow light curves separately, and 
get the respective fitting parameters, such as toy, £xy, and the u> values of each break). 
A minimum number of components (SPL, BPL, or TPL) are introduced based on eye 
inspection of the global features in the lightcurve. If the reduced y 2 is much larger than 
1, we continue to add more components and re-do the fits, until the reduced y 2 becomes 
close to 1 (usually less than 1.5). The reduced y 2 values for some lightcurves are much 
smaller than 1, indicating that some model parameters are poorly constrained. For 
these cases, we fix some parameters and redo the fits until the reduced y 2 becomes 
close to 1. Some GRBs have erratic fluctuations in the lightcurves with small error 
bars, so that the reduced y 2 is much larger than 1. For these cases, we do not add 
additional components to fit the lightcurves, so that their y 2 values remain much larger 
than 1. 

• Next, we jointly fit both optical and X-ray lightcurves by introducing a same tb- We 
search for a possible achromatic break time in the range [toy, ixy]- We still fit the op¬ 
tical and X-ray lightcures at a test break time tb separately in this step. The individual 
y 2 of the optical or X-ray band could not represent the goodness of the jointly fit. To 
evaluate the goodness of the fits for optical and X-ray lightcurves at we introduce 
a weighted reduced y j otal , which is essentially the average reduced y 2 in both bands. 
Taking GRB 050922C as an example: a best join fit is achieved at tb = 17.3 ks, where 
the reduced Xx values are 175/157 and 175/148 ~ 186/157 for the optical and X-ray 
bands, respectively, so that Xtotai can t> e expressed as 361/314. For all the GRBs, we 
search for the common tb with the best Xtotai- We accept the fits with the Xtotai — 3, 
and regard it as not inconsistent with being achromatic^. Usually the parameters of 
this best join fits does not correspond to the best reduced y 2 in each band. 


2 The adoption of a separation line at Xtotai around 3 is somewhat arbitrary, but the value is determined 
based on close inspection of the fitting results of individual bursts. Our results indicate that most GRB 
afterglow light curves are well fit with the BPL or SPL light curves models, with a typical value Xtotai = 1-21± 






















• If both the optical and X-ray lightcurves decay as a SPL, we do not need to search for 
a common break time. The weighted reduced Xtotai is calculated based on the above 
algorithm for the SPL fits in each band. 

• If one band decays as a BPL, while the other band does not have enough data to 
search for a break time and decays as a SPL (e.g., the Grade II or IV in Section 3), we 
impose t b identified in the first band as the common t b , and perform the Xtotai analysis 
as described above. 

The fitted results are presented in Figure Q315J The parameters of the PL or BPL fits 
of all the lightcurves are presented in Table [Q Some lightcurves have additional features 
(e.g., steep decay phase, flares, rebrightening features) in one band. We do not report them 
in Table [0 Our analysis below discards these extra components since they likely arise from 
additional emission components (e.g., in the internal dissipation regions such as internal 
shocks and internal magnetic dissipation sites) other than the external shock. 


3. External Shock Models: Closure Relations and Light Curve Types 

3.1. Closure Relations 

The standard external shock models of GRB afterglows have clear theoretical predic¬ 
tions that can be verified or falsified by the observational data. These models attribute the 
multi-wavelength afterglow emission to synchrotron radiation of electrons accelerated in the 
shock front as the fireball jet interacts with the circumburst medium. The models largely 
do not depend on the details of the central engine activities, so that the afterglow behaviors 
only depend on a limited number of parameters. In the convention of F u oc t -c V -/3 , where a 
and ft are the temporal and spectral indices of the afterglows that can be measured directly 
from observations, the models predict certain relatio nships between a and ft values, which 
are called the “closure relations” of the models (e.g., IZhang fe MeszaroslI2004J: IZhang et al. 
20061 : Gao el al. 20 131). Technically there are many sub-models (e.g., ISM vs. wind, adi¬ 


abatic vs. radiative, whether or not there is energy injection), physical regimes (reverse 


0.50. However, some GRBs (e.g., GRB 050730, 060904B, 080319C, 100901A, 120326A) show a relatively 
large Xtotai > which are around or even slightly larger than 3. Inspecting their light curves, the relatively 
large Xtotai caused by complicated features in the light curves (such as small flares and fluctuations), 
especially in the optical band (e.g., GRB 060904B). However, the PL and BPL fits in any case catch the 
general features of these light curves. Since we are interested in the achromatic/chromatic properties rather 
than the flaring features of the light curves, a relatively loose criterion (Xtotai ~ 3) is reasonable. 
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shock crossing phase, self-similar phase, post-jet-break phase, Newtonian phase), and spec¬ 
tral regimes (different orders among the observed frequency (u) and several characteristic 
frequenci es v r , the self-absorption frequency v a ). We refer to a comprehensive review 


of Gao et ah (120131 ) and references therein for the details of various models. 


For the time frame of our interest (hours to weeks after the trigger), the reverse shock 
crossing phase is usually over, and the blastwave is still in the relativistic phase. This greatly 
reduces the number of relevant models. In Table El we s ummarize the a and f3 predictions 
of various models studied in this paper following Zhang et al.l (hooot) and Gao et al. (2013). 
This includes the ISM and wind models for adiabatic blastwave^], for both pre- and post-jet 
break temporal phases, with and without continuous energy injection, and for two spectral 
regimes (I: v > v c and II: v m < v < u c ) in the slow cooling (v m < u c ) regime. By doing so, 
we have assumed that v a < min (u m , u c ), and min(z/ x , z'o) > which is usually satisfied for 
optical and X-ray afterglow emission for typical GRB parameters. 


2000 


I Him et ah 


The energy injection model invokes either a long-lasting central engine (IDai fe Lulll998l: _ 

Zh ang fe Mesza ros 2001), or a Lorentz-factor-stratified ejecta (IRees fc Meszaroslll998l : ISari fc Meszaros 


2012). The two scenarios are equivalent with each other in terms of 


lightcurve behaviors given a relationship between the central engine parameter q and the 
stratification parameter k (Zhang et ah 2006). We adopt the description of a long-lasting 
central engine with a power-law luminosity history L(t) = Lo(-T)~ g (Zhang & Meszaros 

r ,q 

2001), so that the injected energy is E m j = ■j— 1 t 1 ~ g . The prescription applies when q < 1. 
The relevant closure relations are presented in Table El 

Many observations suggest that GRB outflows are collimated. Assuming a conical jet 
with opening angle 6j, a steepening in the afterglow light curve is predicted when 1/T > 9j 
(T is the bulk Lorentz factor_of the blastwave ). Thu main reason of this steepening is the 


so-called “edge effect” (e.g., Panaitescu et ah 1998)0: The 1/T cone is no longer filled with 


emission beyond the jet break time (when 1/T > 0j). There is a reduction factor in flux 
6]/(1/T) 2 = T 2 9 2 . The relevant closure relations are also presented in Table El 


3 In general, the circumburst medium can be described by an arbitrary profile n oc r~ k . The ISM model 
corresponds to k = 0, and the wind model corresponds to k = 2. In our closure relations, we only consider 
these two cases, since they are naturally expected from the ISM and a pre-explosion stellar wind. For other 
k values, it is not straightforward to imagine a physical mechanism to produce such profiles over a large 
distance scale of interest. We therefore do not include the arbitrary k models in the standard afterglow 
models, but discuss them as possible modified afterglow models. 


4 Sideways ex pansion has been discussed as another factor of steepening the lightcurves ( Rhoads 19991 


Sari et al. 1 9991 . H oweve r, later numerical simulations suggest that this effect is not important (e.g., 


Zhang fc MacFadvenll2009l) . We do not consider this effect in this paper. 
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It is possible that in some GRBs the energy injection phase lasts longer than the jet 
break time, so that a jet break with energy injection both pre- and post-bre ak phases ca n 
be observed. The relevant closure relations of such models were derived in Gao et al. (2013) 
and are also presented in Table © 


3.2. Type of Afterglow Lightcurves 


For the time domain we are interested in and for the optical and X-ray bands, there are 
four types of lightcurves (Fig©: 

(1) Broken power-law lightcurves with an energy injection break: In reference of the 
canonical X-ray light curve (Zhang et al. 2006), as reproduced in Fig©a), the energy 
injection break connects the shallow decay phase (segment II) to the normal decay phase 
(segment III), and a typical light curve is shown in Figure ©b). Before and after the break, 
the adiabatic deceleration a((3 ) relations with and without energy injection (as listed in Table 
[21 IZhang et al.ll2006t iGao et al.ll20131) are used to check whether the data are consistent with 
model predictions. 


(2) Broken power-law lightcurves with a jet break: This corresponds to transition from 
segment III to IV in the canonical lightcurve, and a typical light curve is shown in Figure 
©c) upper curve. Lightcurves of such a category should satisfy the constant-energy, isotropic 
closure relations before the break, and the edge-effect post-jet-break closure relations after 
the break, with no energy injection effect both before and after the break (Table©). The 
post-break decay index is required to be steeper than 1.5 for this model. 


(3) Broken power-law lightcurves with a jet break with energy injection: This model 
allows the energy injection extend to a duration longer than the jet break. The temporal 
break is still defined by the edge effect of a canonical jet, but the decay slopes before and 
after the break are shallower than the previous case (lower curve in Fig©c)), so that a q 
parameter is introduced for both pre- and post-break phases. 


(4) Single power-law decay: For some GRBs, a SPL function is adequate to describe the 
afterglow data (Figure ©d)) after the deceleration phase. In the X-ray band, there might 
be a steeper decay p hase before this SPL phase, which is due to the tail emission from 
the prompt emission (iTagliaferri et al.l 120051; iBarthelmv et al.l l2005t IZhang et al.l 120061 1 . We 
ignore the steep decay phase and treat it as a SPL decay (upper curve of Fig©d)). Similarly, 
in the optical band, some GRBs show an early rising phase, which is a signature of the onset 
of afterglow at the deceleration radius (peak of the lightcurve, lower curve of Fig©d)). We 
treat these lightcurves also as SPL decay ones. 
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For all the types, sometimes there are X-ray flares overlapping on the power-law decay 
segments. We do not include the flares in our data fitting, since they originate from a 


different emission comp onent due to late central engine activities (e.g., IZhang et al.l 12006 


Maxham & Zhang 2009). 


One important task is to perform a self-consistency check between the optical and X-ray 
bands. If a GRB is consistent with the external forward shock model, we demand that the 
GRB satisfies the following criteria: 


• The X-ray and optical lightcurves are consistent with having an achromatic break if 
any; 

• Both the X-ray and optical lightcurves should satisfy closure relations of a same cir- 
cumburst medium type (ISM or wind) in both pre- and post-break temporal segments; 

• Either both bands belong to the same spectral regime, or the two bands are separated 
by a cooling break z/ c , with the X-ray band above the break and the optical band below 
the break (with allowance of a grey zone, see more discussion below); 

• The inferred electron spectral index p from both bands and from both pre- and post¬ 
break segments should be consistent with each other within error; 

• For energy injection models, the energy injection parameter q values derived from the 
X-ray and optical bands should be consistent with each other. 


Technically, we check the consistency between the closure relations for individual tempo¬ 
ral segment in individual energy band. To ensure a same p value derived for different bands, 
we also check the consistency between the data and models in the A/3x,o ~ Aax,o plane. 
Here Aax,o = ccx — cto is the difference between the decay indices in the X-ray and optical 
bands, respectively, in a same temporal segment, and A/3 X: o = fix — fio is the difference be¬ 
tween the spectral indices in the X-ray and optical bands, respectively. Based on the closure 
relations (Table [2]), one can derive the A fix,o ~ Actx.o relations of all the models (Tabled] 
and [4]) . One can see that even though a and fi values can be very different in different mod¬ 
els, the A«x,o and Afix,o values have several well-predicted values. In particular, for the 
SPL, and jet break models, both pre- and post-break values are well-defined constants. For 
the energy injection breaks, the post-break segment does not depend on the free parameter 
q. As a result, if one focuses on the second component only, all the models can be expressed 
as several representative coordinate values in the Afix,o — Aax,o plane. Considering the 
possible grey zones (see below for details), these points define several straight lines in the 
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A fix, o — Actx.o space (FigJTJfor details). If the observed data intersect with these model 
lines (within error), one can regard them as being consistent with the model predictions. 

Take the energy injection break as an example, our analysis uses the following procedure: 
(a) Use the observed spectral indices fio and fix to predict the post-break temporal indices 
«o ,2 and ax ,2 in two possible spectral regimes. Then compare these theoretical predictions 
with the observational values. If theoretical values are consistent with the fitting results 
within error then go to next step. Otherwise, it indicates that this GRB does not fall into 
this light curve type; (b) Use the identified spectral regime to calculate the electron spectral 
index p from the spectral index fi, i.e. p = 2fi + 1 for v m < v < v c , or p = 2 fi for v > v c . 
Compare the p values derived from the optical and X-ray data, respectively. If po = px 
within error, then move to the next step. Otherwise, this GRB does not fall into such a light 
curve type; (c) Use the inferred p value and spectral regimes to infer the energy injection 
parameter q using the temporal index before the break (ao.i and ax,i). Compare the derived 
q values from optical and X-ray bands, respectively. If qo = qx within error, then move to 
the next step. Otherwise, this GRB does not fall into such a light curve type; (d) Using 
the Afix,o — Aax,o relation to double check the data, if the data fall into the predicted 
region in the Afix,o ~ Aa x ,o plane, then this burst can be fully interpreted by such a model. 
Otherwise, the burst does not fall into this category. 


The simplest analytical model (Sari et al. 1998) predicts fi = p/2 for Regime I (y > 
u c ) and fi — (p — l)/2 for Regime II (u m < v < u c ). Detailed numerical calculations 
( Uhm fe Zhang 2014a) showed that the transition between the two regimes may take several 
orders of magnitude in observer time. As a result, some “grey zones”, with (p — l)/2 < 
fi < p/2 are allowed by the model. Therefore the parameter space between the two closure 
relation lines defined by the two spectral regimes in the a — fi plane is allowed by the theory. 
Data points falling into this grey zone should be regarded as consistent with the model. 
There are three possibilities: (1) the optical band is in Regime II, while the X-ray band is in 
the grey zone; (2) the X-ray band is in Regime I, while the optical band is in the grey zone; 
and (3) both bands are in the grey zone. 


For the cases that both the optical and X-ray bands are in the same spectral regime, 
we demand that three spectral indices be the same within error, i.e. fio = fiox = fix , 
where fiox is the spectral index between the optical and X-ray band in the joint spectral 
energy distribution (SED)@. If the two bands are in different spectral regimes, we demand 


5 In order to obtain dox, we roughly fit the SED from optical to X-ray bands. For the optical band, we 
chose the R-band where extinction correction is negligible. For the X-rays, we use the Swift XRT data and 
adopt a typical band 1.5-2 keV, where the absorption effect is negligible. 
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Po < Pox < Px or Pq < Pox < Px- 


4. Confronting Data with Models 

4.1. Grading criteria and sample definitions 

With the above preparation, everything is in place for us to systematically confront the 
broad-band data with the external forward shock afterglow models. Based on how badly the 
data violate the models, we define the following five grades (see also Table [5]): 


• Grade I: Both X-ray and optical bands have SPL lightcurves or BPL lightcurves with 
an acceptable achromatic break. Both bands satisfy closure relations and are self- 
consistent (same medium type, p and q values). These are the best examples where 
the GRB afterglow data abide by the external shock model predictions; 

• Grade II: Some GRBs have a clear break at t h in one band (e.g., X-rays), but do 
not have a break in another band (e.g., optical). The missing break is likely due to 
incomplete observational coverage before or after the break. The data are consistent 
with the hypothesis of an achromatic break, and both bands satisfy closure relations 
self-consistently. These GRBs are almost as good as Grade I in terms of abiding by 
the external shock models; 

• Grade III: Both X-ray and optical bands have SPL lightcurves or BPL lightcurves with 
an acceptable achromatic break. However, at least one temporal segment in one band 
does not satisfy the closure relations in a self-consistent manner with respect to other 
segments/band. 

• Grade IV: This is the Grade II equivalent for Grade III. One band does not have a 
break, but the data are consistent with the hypothesis of having an achromatic break. 
At least one temporal segment in one band does not satisfy the closure relations in a 
self-consistent manner with respect to other segments/band. 

• Grade V: Clear evidence of chromatic breaks and violation of closure relations. These 
GRBs cannot be interpreted within the one-component external shock model^l 


6 Some of these GRBs may be still in terpreted within two- component external shock models with each 
component dominating one band (e.g., De Pasauale et al. 20091 1. However, the demanded parameters for the 


two components are rather contrived. 
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With these five grades, we define three samples: 


• Gold Sample : The GRBs in Grade I and II are defined as the Gold sample GRBs, since 
no observed information violates any predictions of the external shock models; 


• Silver Sample : The GRBs in Grade III and IV are included in this sample. Even though 
at least one segment/band does not satisfy the closure relations self-consistently, the 
basic requirement of achromaticity is not violated. We note that the closure relations 


plicated models invoking, e.g., a structured jet ( 

Zhane & Meszaros 

2002; 

Rossi et al. 

2002; 

Kumar & Granot 2003 

; Granot & Kumar 2003) or a circumbur 

st density medium 
redict light curve 
urthermore, if the 

with an arbitrary k value (at least for a certain distance range), p 
behaviors that may not fully abide by the simple closure relations. F 


GRB engine is long-lived and a long-lasting reverse shock outshines the forward shock, 
a variety of rich light c urve behaviors can be generated, which do not follow the simple 


closure relations (e.g., lUhm et al .1120121; lllhm fc Zhandl2014bl ). So it is possible that 


the GRBs in the silver sample are still consistent with the external shock models; 


• Bad Sample: The GRBs in Grade V violate the basic achromaticity principle of the 
external shock models and do not abide by the closure relations, and therefore cannot 
be interpreted within the framework of the external shock models. 


4.2. Grading results 

The 85 well-sampled GRBs in our sample are graded based on the above-defined grading 
criteria. The GRBs in the five grades are presented in Figuresdj-O respectively. The relevant 
data of different grades are presented in Table Q] and [6] 

• Grade I: As can be seen from Table [Q and [6J. and Figure [lj within errors 43/85 GRBs 
satisfy the Grade I criteria. Out of 43 GRBs, 13, 8 and 22 GRBs are constrained to 
have an energy injection break, jet break and SPL decay, respectively. 

• Grade II: within error 2/85 GRBs fall into this grade (Figj2|). 

• Grade III: there are 34/85 GRBs falling into this grade (Figj3j). Among the sample, 
15/34 and 19/34 GRBs have SPL and BPL lightcurves, respectively. GRBs 060906, 
080319B and 100219A have two beaks at different times, respectively. 

• Grade IV: there are 3/85 GRBs falling into this grade (Fig]4j). 
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• Grade V: there are 3/85 GRBs falling into this grade (Figj5]). Two of them (GRBs 
060607A and 070208) show clear chromatic breaks with good temporal coverage in 
both bands at the break times. One GRB (GRB 070420) shows a chromatic behavior 
based on the available data and simple model fitting, even though no observational 
data are available in the optical band at the break time of the X-ray band, so that the 
existence of a break in the optical band (even though very contrived in shape) at the 
same epoch is not completely ruled out. 

Consequently, we get three samples: 

• Gold sample: This sample has 45/85 GRBs, including 13/49, 8/49 and 24/49 GRBs 
satisfying the energy injection, jet break, jet break with energy injection, and SPL 
decay models, respectively. Among them, 27/49 and 18/49 are consistent with the 
ISM and wind models, respectively; 17/49, 4/49 and 24/49 GRBs are consistent with 
being in a same spectral regime, different spectral regimes (X-ray band in regime I and 
optical band in regime II), and grey zone, respectively. Among the 17 GRBs with the 
same spectral regime, 15 and 2 GRBs are consistent with being in the ISM II and wind 
II spectral regimes, respectively. For the 4 GRBs with different spectral regime, all of 
them are consistent with having an ISM medium. 

• Silver sample: This sample has 37/85 GRBs, which may (or may not) be interpreted 
within the more complicated numerical external shock models. 

• Bad sample: Only 3/85 GRBs definitely violate the basic achromaticity principle of 
the external shock models and therefore belong to the bad sample. 

Figure [7^ shows A/?x,o — Aax,o distributions for Gold sample. For the energy injection 
sample, we only used the post-break segment to remove the g-dependence. These are the 
GRBs that also satisfy the closure relations in all temporal segments. We do not show the 
closure relation a — (3 plots since the energy injection models have an extra g-dependence on 
the a values. To show the details of how each burst may fall into the model predictions of 
each model (grey zone included), in Figure [7)(b-e) we show the A/3x,o ~ Aax,o distributions 
of those Gold-Sample GRBs that satisfy the ISM and wind medium models with p > 2 and 
1 < p < 2, respectively. The Silver sample GRBs are collected in Figure [7f). About half of 
them fall outside the predicted region (red box) defined by the models. Even though some 
fall into the box, they do not satisfy the closure relations in all the temporal segments in all 
energy bands. 


5. Statistics of the External Shock Afterglow Model Parameters 


Since the Gold sample (Grade I and II) GRBs comply with the external shock models 


well, they serve as an excellent sample to study external shock model parameters. The 
derived external shock parameters of the Gold sample GRBs are presented in Table [6] We 
present some statistical properties of these model parameters in this section. 


5.1. Temporal indices a 


Figure [8] shows the distributions of the temporal indices a in different energy bands 
and different temporal segments. They are all well fitted with Gaussian distributions for 
each band/temporal segment. For the GRBs having a BPL lightcurve, the typical a values 
are «o,i = 0.49 ± 0.45, op ,2 = 1-44 ± 0.39, ax, 1 = 0.58 ± 0.63 and ax ,2 = 1-50 ± 0.27, 
respectively (FigJHK)• For the GRBs with a SPL lightcurve, one has ao = 1.26 ± 0.38, and 
ax = 1.39 ±0.26 (Figj8t>). For the BPL sample, we also separate it into the energy injection 
sample and the jet break sample and perform the statistics. For the energy injection breaks, 
one has «o,i = 0.25 ± 0.12, op ,2 = 1-26 ± 0.26, ax, 1 = 0.30 ± 0.27 and ax ,2 = 1-35 ± 0.24, 
respectively (FigJSfc). For the jet breaks, one has «o,i = 0.77 ± 0.18, op ,2 = 1-66 ± 0.16, 
ax, 1 = 0.95 ± 0.16 and ax ,2 = 1-70 ± 0.19, respectively (FigJHJi) . Both the pre-break and 
the post-break a values in the energy injection sample are systematically shallower than 
those in the jet break sample. On average, the X-ray lightcurves are steeper than the optical 
lightcurves, consistent with the expectations of the theoretical models (i.e. the X-ray band 
is more likely above v c while the optical band is more likely below u c ). 

Another self-consistency check is to compare the observed change of decay slope, A a = 
«2 — «i, with the model predictions. From the closure relations (Table |2j), one can derive 

• For energy injection breaks: 


( (l-g)(2+/?) 


ISM II (p > 2) 

ISM II (1 < p < 2) , 

ISM I (p > 2), wind II (p > 2), wind I (p > 2) 


2 ’ 


(i—g)(i9+2/3) 



Aa ~ < (1 — qfa+P) ’ 


ISM I (1 <p< 2) 
wind II (1 < p < 2) 
wind I (1 < p < 2) 


(5) 



4 


(l-q)p3+/9) ’ 

4 ’ 


• For jet breaks: 
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f |, ISM I and II (p > 2 and 1 < p < 2) 
\ -f, wind I and II (p > 2 and 1 < p < 2) 


• For jet breaks with energy injection: 


( 6 ) 


ISM II (p > 2), ISM I (p > 2 and 1 <p < 2) 

ISM II (1 < p < 2), (7) 

l, | , wind II and wind I (p > 2 and 1 < p < 2) 

Figure [9] shows a comparison between the observed Aa 0 bs and the theoretically predicted 
Aa t h for each GRB derived from the measured (3 and q values using the corresponding closure 
relations. One can see that the two are consistent with each other. 

Figure [TO] displays the observed Act distributions of various samples. For the Gold 
sample, the Act distributions of optical and X-ray data are consistent with each other, i.e. 
Acto = 0.94 ± 0.23 and Actx = 0.88 ± 0.28 fFigfTOhh Furthermore, the Act values of both 
bands in sub-groups (energy injection breaks and jet breaks) are also consistent with each 
other: Acto — 1-05 ± 0.17 and Actx = 1-10 ± 0.21 for the energy injection breaks, and 
Acto = 0.75 ± 0.22 and Actx = 0.75 ± 0.22 for the jet breaks (FigITOb). The Silver sample, 
on the other hand, shows a poorer statistical behavior fFigJTUb and d). 


5.2. Spectral indices (3 

Figure |TT] shows the spectral index distributions for the Gold sample. In general, the 
distributions can be fitted with gaussian functions. For the global sample, one has (3q = 
0.70 ± 0.15, and f3x = 0.98 ± 0.15 fFigfTThh In the Gold sample, 17/45 GRBs have both 
the optical and X-ray bands in the same spectral regime. One has j3o = 0.77 ± 0.19, and 
/ 3x = 0.89 ± 0.15, which are consistent with each other (FigJTTJo). The rest 28/45 GRBs are 
identified to have X-ray and optical bands separated by a cooling break. The results show 
flo = 0.68 ± 0.18, (3x = 1-01 ± 0.14, with A (3 = (3x — (3o = 0.37 ± 0.18, which is consistent 
with the theoretically expected value 0 < A (3 < 0.5 (FiglTTb). 

We investigate the /3 distributions in different types of lightcurves. For the energy 
injection sample, one has (3o = 0.78 ±0.12, and (3x = 1.01 ±0.13 (FiglTTld): for the jet break 
sample, one has (3q = 0.59 ± 0.11, and /3x = 0.97 ± 0.08 fFigJTTbk and for the SPL sample, 
one has f3o = 0.74 ± 0.24, and /3 X = 0.95 ± 0.19 (TigfTTlf h 

We also investigate the (3 distributions in different ambient medium types. For the ISM 
model (27/45 GRBs), one has /3q = 0.72 ± 0.21, and f3 x = 0.98 ± 0.10 fFigfTTkk and for 
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the wind model (18/45 GRBs), one has (3q = 0.70 ± 0.10, and /3x = 1-00 ± 0.20 (FigJITli). 


results (e.g., 


Panaitescu & Kumar 

2002; 

Yost et al. 

p _ '' j 

2003; 

Zhang et al. 

2006; 

Schulze et al. 


2011 ). 


5.3. Electron spectral index p 


Figure H2l shows the distributions of the electron spectral index p of the Gold sample. 
It has a Gaussian distribution with p = 2.33 ± 0.48 (FigJT2h). which is very consistent 
with the typical value of p for relativistic shoc ks due to lst-order Fermi acceleration (e.g., 


Achterberg et al. 

2001; 

Ellison & Double 

2002) 

It also has a wide distribution, which is 

consistent with previous studies (e.g., 

Shen et al 

2006; 

Liang et al. 

2007, 

2008; 

Curran et al. 


20101 ). 


The p distribution in different sub-samples are also generally consistent with each other. 
Within the Gold sample, those GRBs with optical and X-ray bands in the same spectral 
regime have p = 2.58 ±0.39, whereas those with optical and X-ray bands in different spectral 
regimes have p = 2.17 ± 0.44 (FigJl2h). For the three light curve sub-samples, one has 
p = 2.34 ± 0.38 for the energy injection break sample, p = 1.91 ± 0.37 for the jet break 
sample, and p = 2.48 ±0.47 for the SPL sample, respectively (Fig]T2b). For the two medium 
type models, one has p = 2.43 ± 0.57 for the ISM model, and p = 2.28 ± 0.33 for the wind 
model, respectively (FigJT2fc). 


5.4. Break time R 

Figure [IS] shows the distributions of the observed achromatic break times, R. The global 
distribution in the Gold sample gives log(R/ks) = (3.8±0.9). Separating the energy injection 
sample and jet break sample, one has log(R/ks) = (3.6 ± 1.9) for the energy injection break 
sample, log(R/ks) = (3.9 ± 0.7) for the jet break sample. The energy injection ends (which 
depends on central engine) is on average earlier than the jet break time (which depends on 
geometry of the jet). The distribution of the energy injection break time is wider than the 
jet break time distribution. 
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5.5. Energy injection parameter q 

Within the Gold sample, 13/45 GRBs show an energy injection type break, i.e. either 
an energy injection break or a jet break with energy injection. Among them, 4/13 and 
9/13 GRBs satisfy the ISM and wind model, respectively. The distributions of the energy 
injection parameter q of various samples are shown in Figure dH The global sample has 
q = 0.22 ± 0.11. The ISM and wind models have q = 0.20 ± 0.12 and q = 0.23 ± 0.13, 
respectively, which are consistent with each other. 


5.6. Shock magnetic field equipartition factor eg 


Among the derived shock parameters, the magnetic field equipartition factor is of 
special interest. If the shock simply compresses the upstream magnetic field, then the ex¬ 
pected 6b is low, of the order of 10~ 6 — 10~ 7 . If, however, various plasma instabilities are 


playi ng a role to amplify the magnetic fields (e.g., iMedvedev fc Loeblll999l : iNishikawa et al. 


2009 ), one would expec t a relatively large 6 b as high as 0. 1 . Early afterglow modeling (e.g., 


Wiiers fe Galamal Il999l: IPanaitescu fe Kumarl l200ll . 12002c lYost et al.l 1200311 derived a rela¬ 
tively large 6b, with a typical value ~ 0.01. On the other hand, modeling of GeV emission 
in several Fermi/LAT-d etected GRBs led to the suggestion t h at 6b should be relati vely low 
at least for some GRBs (IKumar fe Barniol Duranll2009l. 1201011 . ISantana et al.l (1201411 derived 
6b for a large sample of GRBs, and derived a medium value of 10 -5 . 


With our Gold sample, we can constrain 6b independently. Even though in most GRBs, 
6b cannot be constrained due to the degeneracy of the data, one can still place interesting 
upper limits to 6b based on the medium type and spectral regime of the GRBs. For example, 
in the ISM model v c decreases with time. For a regime II (u m < v < z/ c ) GRB, the last data 
point in the light curve would set a lower limit on v c at that epoch, and hence, an upper 
limit on 6b- Similarly, in the wind model v c increases with time. For a regime II GRB, the 
first data point in the light curve would set a lower limit on v c at that point, and hence, an 
upper limit on 6b- 


In the Appendix, we present expressions of u m , u c , max and the kinetic energy of the 
afterglow, Ey isn . For the v > 2 cases, we ad opt the formalism in previous works (IZhang et al. 


2007al; iGao et al.l 120131: iLii fe Zhang I2014T) . New expressions for the 1 < p < 2 regime are 


also presented following [Gao et ah (12013n . We then derive the expressions of 6b in various 
models, Eqs. (IA6D . (IA9D . (IA20D . and (IA23D . which are used to constrain e#. 


The derived upper limits of 6b are presented in Figure [15j with other parameters fixed 
as 6 e = 0.1, n = 1 or A* = 1. One can see that in general these upper limits point towards 
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a relatively low 6b value. In some cases for the ISM model, the upper limit s are eve n lower 
than 10 -5 . These resu lts ar e cons istent with the finding s of lKumar fc Barniol Duranl (12009 


2010 1. ISantana et al.l ( 2014T) . and lBarniol Duranl ( 20141) . 


5.7. Energetics 

The isotropic y-ray energy E 7jiso is calculated as 


tp _ 4 nDlS^k 

7 ’ iso 1 + z ’ 


( 8 ) 


where S 1 is the gamma-ray fluence in the BAT band, D l is the luminosity distance of the 
source at redshift z, and the parameter k is a factor to correct the observed y-ray energy in 
a given band pass to a^ broad band (e.g., 1 — 10 4 keV in the rest frame) with the observed 
GRB spectra (IBloom et ahl 20 011 ). It is well known th at a typical GRB spectrum is well 
fitted with the so-called Band function ([Band et ah 1993 ). If the Band parameters are 
measured for a burst, these parameters are used to derive the k parameter. However, owing 
to the narrowness of the Swift/BAT band, the spectra of many Swift GRBs in our sample 
are adequately fitted with a single power-law, N oc E ~ l , so that the Band parameters 
are not well constrained. F or these GRBs, we use an empirical relation between E r and 


BAT-band photon index T ( Zhang et ahl l2007bl: ISa.ka.moto et ahl 120091: IVirgili et ahl 1201211 
to estimate E v . Takin g typ ical values of the photon indices a = —1.1 and /3 = —2.2 


(IPreece et al.ll2000t iKaneko et al.ll20061 ). we can derive E 7 ; so values of the GRBs with redshift 
measurements in our Gold sample, which range from 10 51 to 10 55 erg, with a typical value 
log(E 7jiso /erg) = 53.15 ± 0.69 (FigH^a)). 

The isotropic kinetic energy of the afterglow E K [ so can be derived from the afterglow 
data. In gen eral, broad-band modeling is needed to precisely measure EK,iso (jPanaitescu & Kumar 


2001, 200211 . Most GRBs do not have adequate data to perform such an analysis. More 
conveniently, one may use the X-ray data only to constrain E K iso , since the X-ray band 


only weakly depends on e e ( 

Kumar 

1999; 

Freedman & Waxman 

2001; 

Berger et al. 

2003; 

Lloyd-Ronning & Zhang 

2004; 

Zhang et al. 

2007a 

), see Appendix for detailed derivations. 


For the cases with energy injection, E -&,iso is a. function of time during the energy injection 
phase. Following [Zhang et ahl (I2007al f. we calculate E K iso at two different epochs, one at 
the break time 4, when energy injection is over, and another at a putative deceleration 
time tdec ~ max(60 s, T 90 ). We use the X-ray flux at U to derive EK,end, and then derive 
Ek, dec = -EK.end^dec/G) 1-9 - The total injected energy is calculated as E K}inj = E K , e nd-^K,dec- 
Based on our constraints on 6b, we take 6b = 10~ 5 for all the GRBs in our calculations. 
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Other parameters are taken as typical values: e e = 0.1, n = 1 or A * = 1, and Y = 1. 

In Figure [TBl we present several statistical results of the K K ,i S o calculations. For the 
GRBs with energy injection, the distributions are log (Ek, end/erg) = 54.99 ±0.86 (FigJIBb). 
log (±K,inj/erg) = 54.95 ± 0.61 (FigJl 6 b). and log (K k , dec/erg) = 53.29 ± 0.45 (FigJl 6 H). For 
the entire Gold sample, one has log (Ek, dec/erg) = 54.66 ±1.18 (FigJl 6 ld). It is interesting to 
see that the energetics of the energy-injection sample reache a similar level as the no-energy- 
injection sample after the energy injection is over. Clear correlations are found among 
different energy components: -Effiend - hk.inj relation E K ,inj ,52 = 0.69E^° 2 d 0 5 2 2 (FigUZk), 
Ek,dec Ek ,inj relation /.k. ini . n 2 41.7Ek j ec ^ (FigUTjo), E 7i i SO -Ek,end relation Expend ,52 

476.3^f s ^2 23 (FigUTfc). E 7i i SO — EK,dec relations in the entire Gold sample, ±K.dec,52 = 
56.2S 7 0 J 3 o f 5 ° 2 19 (FigUZH), in the energy injection sample, E K , d ec ,52 = 8.9E}£3£ 29 (FigEb), 
and in the no-energy-injection sample, ±K,dec ,52 = 316.2E^f s ^2 21 (FigJTTf). respectively. In 
particular, the±K,inj ; 52 = 41.7^’^52° correlation suggests a substantial energy injection 
during the shallow decay phase for most GRBs. 


5.8. Radiative efficiency ?] 7 


The GRB radiative efficiency, defined as (Lloyd-Ronning & Zhang 200411 


ih 


v 7,iso 


-E'7,iso “ 1 “ -E'K,iso 


(9) 


is an essential parameter to probe how efficient a burst converts its global energy to prompt 7 - 
ray emission. As mentioned above, if there is continuous energy injection, the kinetic energy 
of the afterglow Ek,i SO takes different values if one chooses different epochs. In principle, ?? 7 
can be defined for two different e poc hs, td ec and t&, which have different physical meanings 
(see a discussion in Zhang et al. 2007a). 


Figure fl8l shows the radiative efficiencies calculated at fdec and G as a function of E 7iiso 
along with their histograms. No significant correlation between rjy and E 7 i i SO is found. The 
fireball internal shock model predic t s a relatively small effic i ency of a few per cen t (e.g., 

Maxham fe Zhang 2009 : Gao fe Meszaros 20151 b Pre¬ 


Kumar 

1999; 

Panaitescu et al. 

1999; 


vious cons traints on GRB radiative efficiencies give relatively large values, as large as above 


90% (e.g., iLlovd-Ronning fc Zhanel 12004 IZhang et al.l l2007al: iRacusin et al.l l 201 ll ) for some 


GRBs. This challenges the internal shock models, and favors alternative prompt emission 
models, such as diss ipation of magnetic fields (IZhang fe Yanll201ll) or photospheric emission 
( Lazzati et ah 2013 k 


The derived efficiencies can be fit with rough log-normal distributions fFigJTBb. d, f). For 
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the entire sample, one has log(7/ 7 i dec/%) = 0.75 ± 0.86. For the sub-sample GRBs without 
energy injection, the radiative efficiency is lower, with log(?/ 7 i dec/%) = 0.37 ± 0.61. For 
the sub-sample GRBs with energy injection, the radiative efficiencies read log(r] 1:dec /%) = 
0.92 ± 1.25 for t dec , and log(77 7iend /%) = —0.89 ± 0.97 for t bl respectively. 


The deri ved efficiencies are somewhat smaller than the values derived in previous work 
(e.g. [Zhang et ah 2007 a). The main reason is the adoption of a smaller value of e# ~ 10~ 5 , 
so that the derived E^^ iso are systematically larger. This greatly alleviates the low-efficiency 
problem of the internal shock models. Nonetheless, some GRBs still have te ns o f percent 
efficiency, which demands a co ntrived setup for the internal shock models fe.g.. iBeloborodov 
2000l; iKobavashi &; Saril l200ll ) . If fdec is a dopted, which is m ore natural for most prompt 


emission model to calculate efficiency (see IZhang et ah 2007a for a detailed discussion), ?/ 7 


is still typically too large for the internal shock model. This is on the other hand consistent 
with the suggestion that internal collision-induced magnetic reconnection and turbulence 
(ICMART) is the dominant process to power GRB prompt emission in the majority of 


GRBs, whi ch can typically gives tens of percent radiative efficiency flZhang fe Yanl 2011 


Deng et al.l I2015T ) . This conclusion is also consistent with i ndependent studies of modeling 


the GRB prompt emi ssion spectrum flllhm fe Zhang|l2014d ) and quasi-thermal photosphere 
emission component f Gao &; Zhane 2015h . 


5.9. Jet opening angle and geometrically-corrected gamma-ray energy 


In the Gold sample, 8/45 GRBs show a jet break. These include five GRBs of without 
energy injection and one more with energy injection. The ambient medium type of all 6 
GRBs is ISM. Five out of these six GRBs have redshift information (Table [7]). 


Under the ass umption of a conical jet, one can derive the jet opening angle based on 
observational data (IRhoadsIIl999l: ISari et al.lll999l: iFrail et al.ll200ll) : 


0,- = 0.070 rad 


h 


d 


3/8 


1 day/ 


1 + z 


-3/8 


E] 


K,iso 


IQ53 er g s 


- 1/8 


n 


0.1 cnr 3 


1/8 


( 10 ) 


We then calculate the geometrically corrected 7 -ray energy 


Ery ( 1 COS 0 j ) 5 


( 11 ) 


and kinetic energy 


E K = (1 - COS0j)£K,iso- (12) 

Here -E K ,i S o is taken as E Kend for the energy injection sample. The medium density is taken 
as n = 1 cm -3 .The results are presented in Table [7] and Figure [191 The best fitting results 
give 9j = (2.8 ± 1.5)°, log(U 7 /erg) = 49.86 ± 0.65, and log(U K /erg) = 50.89 ± 0.54. 











































6. Conclusions and Discussion 


The chromatic afterglow behavior observed in some GRBs has raised the concern regard¬ 
ing whether the external forward shock models are still adequate to interpret the broad-band 
afterglows of GRBs and whether alternative ideas, e.g., a long-lasting engine-driven after¬ 
glow, are needed to account for the data. In order to answer “how bad/good the external 
shock models are”, in this paper, we systematically studied 85 Swift GRBs up to March 
2014, which all have high-quality X-ray and optical light curves and spectral data to allow 
us to study the compliance of the data to the external forward shock models. The results of 
this study can be summarized as the following. 

Based on how well the data abide by the external forward shock afterglow models, we 
categorized GRBs into five grades and three samples: 

• A Gold sample (Grade I and II) includes 45/85 GRBs. These GRBs are fully consistent 
with the theoretical predictions of the external shock models, including having an ac¬ 
ceptable achromatic break and fulfilling various closure relations between the temporal 
decay indices a and spectral indices f3. 

• A Silver sample (Grade III and IV) includes 37/85 GRBs. These GRBs are also consis¬ 
tent with having an acceptable achromatic break, even though one or more afterglow 
segments do not comply with the closure relations. These GRBs are potentially inter¬ 
pretable within the framework of external shock models. 

• A Bad sample (Grade V) only includes 3/85 GRBs. These GRBs show direct evi¬ 
dence of chromatic behaviors, which cannot be accounted for within single-component 
external shock models. 

The bottom line of this study is to address how bad/good the external shock models 
are. Our results show that external shock models work very well for at least ~ 53% of GRBs 
(our Gold sample). These GRBs can be interpreted within the simplest afterglow models. If 
more advanced modeling invoking other factors (e.g., structured jet or long-lasting reverse 
shock) is carried out, up to ~ 96% of GRBs (including the Silver sample) may be accounted 
for within the external shock models. Only less than 4% GRBs truly violate the basic 
expectations of the external shock models, and demand another emission component (e.g., 
central engine afterglow) to account for emission in at least one band (e.g., the X-ray band). 

Several caveats deserve mentioning. First, we only focused on the main afterglow com¬ 
ponents (SPL, BPL or TPL) of the X-ray and optical lightcurves. In some GRBs, there are 
additional components overlapping with these main components, such as the X-ray steep 
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decay phase, X-ray flares, and optical re-brightening features, which are not included in 
the analysis. These features are usually chromatic, and demand additional emission com¬ 
ponents to interpret the data. The general conclusion that a lot of GRBs have extended 
central engine activities (Zhang et al. 2006) remain valid. The true duration of the GRB 
central engine activities may be much longer than what is measured by the GRB duration 
T 90 flZhang et all 2014). Second, we adopted a relatively “loose” criterion ( Xtotai — 3) to 
define “achromaticity” by requiring the X-ray and optical lightcurves to have a same break 
time. Searching for break times independently in the two bands often results in somewhat 
different break times, but many GRBs can be made being consistent with achromatic. The 
relatively large Xtotai ~ 3 in some GRBs is mostly caused by the additional features (small 
flares and fluctuations which we do not care) in the otherwise (broken) power law lightcurves. 
We therefore believe that our approach is appropriate to address the question of “how bad 
the models are”. On the other hand, if in the future high-quality data indeed show slight 
chromatic behaviors with high confidence, one should take cautious to the fraction numbers 
presented in this paper, and consider how such slight chromatic behaviors may impact the 
models. Finally, we only studied 85 GRBs that have both bright X-ray and bright optical 
emission data to allow us to perform the test. There are more GRBs detected by Swift 
(~ 900 with X-ray lightcurves and ~ 260 with optical lightcurves). Due to the complicated 
sample selection effects, we do not guarantee that the fractions of Gold, Silver and Bad sam¬ 
ples are reliable numbers for the entire GRB population. In any case, 85 GRBs represent a 
reasonably large sample, so that our statistics are valid at least for the “bright” sample of 
GRBs. 


With the Gold sample, we further performed a series of statistical analyses of various 
observational properties and model parameters. Following interesting conclusions can be 
drawn: 


• Temporal index or. The temporal indices a in different bands and different temporal 
segments satisfy the afterglow model predictions. On average, the X-ray lightcurves are 
steeper than optical. For BPL lightcurves, the degrees of the break, A a, are consistent 
with the theoretical predictions of the energy injection models or jet break models; 

• Spectral index j3\ The spectral indices /3 in the optical and X-ray bands are f5 0 = 
0.70±0.15, (3x = 0.98±0.15, respectively. Some (17/45) have X-ray and optical bands 
in the same spectral segment, while most (28/45) have the two bands separated by v c 
or in the grey zone. Statistically, A (3 = 0.37 ± 0.18 is consistent with 
value 0-0.5, a range of A/3, including those expected in the grey zone 
2014a). 

• Electron spectral index p\ The typical value p = 2.33 ±0.48 is very consistent with the 


the theoretical 
f Uhrri fe Zhang 
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theoretical predictions for relativistic shocks. A wide range of p values are observed, 
which is consistent with previous findings. 

• Break time The typical break time is found to be log(G/ks) = 3.8 ± 0.9. The break 
time of energy injection sample (log(i&/ks) = 3.96 ± 1.9) is statistically earlier than 
that of the jet break break sample (log(4/ks) = 3.9 ± 0.7). 

• Energy injection parameter q: the central value is q = 0.22 ± 0.11, and the ISM and 
wind models are consistent with each other, with q = 0.20 ± 0.12 and q = 0.23 ± 0.13, 
respectively. 


• Magnetization parameter es- The derived upper limits of es suggests that the typical 
value of this parameter is low (say, 10~ 5 ), which is consistent with previou s work 
( Kumar &; Barniol Duran 2009 . 2010 : Santana et al. 2014 : Barniol Duran 2014 1. 


• Energetics: The typical isotropic y-ray energy is log(F 7iiso /erg) = 53.15 ± 0.69. For 
the energy injection case, the typical isotropic kinetic energy in the blastwave is log 
(Ek, dec/erg) = 53.29±0.45 at the deceleration time, and log (i K: end/ er g) = 54.99i0.86 
when energy injection is over. For GRBs without energy injection, the typical blastwave 
kinetic energy is log (Ek, dec/erg) = 54.66 i 1.18. Clear correlations among various 
energy components are found. 

• Radiative efficiency y 7 : With a small e# ~ 10 -5 adopted, the derived radiative effi¬ 
ciency ?? 7 is lower than previous studies. For the entire Gold sample, log(y 7i d ec /%) = 
0.75 i 0.86. Yet, the efficiency is still large for some GRBs, especially the ones with 
energy injection. For these GRBs, the efficiency measure at the deceleration time has 
log(7? 7) dec/%) = 0.92 ± 1.25, which still challenges the internal shock model. 


• Jet opening angle 6y. For the jet break sample, we derived the typical jet opening 
angle as 0 3 = (2.5 ± 1.5)°. The jet-corrected y-ray energy and kinetic energy are 
log(F 7 /erg) = 49.86 ± 0.65 and log(F K /erg) = 50.89 ± 0.54, respectively. 
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A. Expressions of Ek,i S o and e# 

In this Appendix, we present expressions of EK,iso and es of the external forward shock 
afterglow models. 


A.l. The ISM model 


For the p > 2 c ase, the forward shock emission can be characterized as (lYost et al.ll2003 
Zhang et al. 20078 ) 


U m = 3.3 X 10 12 HZ (--- ) (1 + ^) 1/2 4 /2 - 2 e e,-l^K , !L,52^ 3/2 5 


p — 1 


Vc = 6,3 x 10 15 Hz (1 + + Y)-h- B 3 p 2 Ep£ i2 n-Hp 12 , 


= 1.6 mjy (1 + z)D 2 p i p_ 2 E K j a , fi2 n 1, ' J , 


d/2 


(Al) 


where E Ki i SO ,52 is the isotropic kinetic energy (in units of 10 52 erg) in the blastwave, 0 is 
the time since trigger (in units of days), n is the density of the constant ambient medium, 
Dl = 10 28 cm D 2 s is the luminosity distance, and 


Y = [-1 + (1 + ^ViVzCe/ £ b ) 1/2 ]/2 


(A2) 


is the inverse Compton parameter, with rji = min[l, (r'c/r'm)^ 2 p ^ 2 ] (Sari & Esin 2001), and 
r /2 < 1 is a correction factor introduced by the Klein-Nishina correction. 


For p > 2 and in the u > max(i/ m , u c ) regime, one has OZhang et al.ll2007aT) 


uF v (u = 10 18 Hz) = F Vt 


l/2 I/ (p-l)/2 z/ (2-p)/2 


maxb On 


X 


= 5.2 x 10 -14 erg s -1 cnF 2 D^( 1 + z) {p+2)/i { 1 + Y)~ x f pl e# /4 e p ~\ 

p(p+2)/4 (2-3p)/4 (2-p)/2 / A o\ 

X A K ,iso,52 t d Os i (Ad) 

where vf„{v = 10 18 i7z) is the energy flux at 10 18 Hz (in units of ergs s -1 cm -1 ), and 

p -1 


f P i = 6.73 


p — 2 
p — 1 


(3.3 x io -6 )^ -2 ' 3)/2 


(A4) 


is a function of electron spectral index p (Zhang et ah 2007a). One can then derive 
uFju = 10 18 Hz) x 4/(p+2) 


v K,iso,52 — 


5.2 x 10 -14 erg s -1 cm -2 

4(l-p)/(p+2) ,(3p-2)/(p+2) 2(p-2)/(p+2) 

Xfc e,-1 L d U 18 i 


A> 28 (P+2) (1 + -) _1 (1 + Y) 4 / (p+2) / p ) 4/(p+2) e^“ p 2 )/(p+2) 

(A5) 
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CB,-2 — 


6.3 x 10 15 Hz 


Vr 


(p+2)/(p+4) 


/zF,(zy = 10 18 Hz) 


-2/(p+4) 


D 


-4/(p+4) 

28 


5.2 x 10 _14 erg s _1 cm _2 / 

X (1 + yj-2(p+3)/(p+4) n -(p+2)/(p+4)y 2 /(p+4) e 2(p-l)/(p+4)^2p/(p+4)^(2-p)/(p+4) 

For p > 2 and in the u m < u < u c regime, one has 

— 10 18 Hz) = F^^-^Vi 3 -^ 2 

= 6.5 x 10- 13 erg s " 1 cm " 2 D£(l + 

xn 1 / 2 4 3 “ 3 p)/ 4 ^ 8 _p)/2 , (A 7 ) 


so that 


A'K,iso,52 


uF v {u = 10 18 Hz) 


4/(p+3) 


t-)8/(p+3)/i . \ — 1 /•—4/(p+3) — (p+l)/(p+3) 

-^28 G "r G J pi t B,~ 2 


6.5 x 10 _13 erg s _1 cnr 2 / 28 v ' y Jpl 

v ,4(1—p)/(p+3) —2/(p+3)y.(3p—3)/(p+3) 2(p-3)/(p+3) 

Afc e,-1 n L d ,y 18 ) 


(AS) 


hB,~2 — 


6.3 x 10 i5 Hz 


z/ r 


(p+3)/(p+4) 


z/F^zz = 10 18 Hz) 


-2/(p+4) 


D 


—4/ (p+4) 


28 


6.5 x 10 -13 erg s _1 cm -2 

y ^ _|_ y^-2(p+3)/(p+4)^ —(p+2)/(p+4)y2/(p+4)^2(p-l)/(p+4)^-2p/(p+4)^(3-p)/(p+4) 


Following Gao et al. ( 20131 ), below we derive the expressions in the p < 2 case. 

In the v > max(p m , v c ) regime, one has 

F K ,iso = {vF v {v = 10 18 Hz)) 16/(p+14) F 32/(p+14) /p 2 16/ ( P+14) n ( P“2)/(P+14) t (3p+10)/(p+14)^ 1 + y)16/(p+14) 

x(l + (Al0) 


f p2 = 0.00529 x e°' 767p 


V 


p 


(All) 


e B = 9.44 x 10 27 z/- 2 / 3 (z/F,(z/ = 10 18 H z))- 16 / 3 (P + 14 )D- 32 / 3 ( P + 14 d^ / 3 (p+ 14 ) n-( 3 P + 26 )/( 3 P +42 ) 

xt -4(p+6)/3(p+14)^ + y^-4(p+18)/3(p+14)^ + ^-(p+2)/3(p+14) e 16/3(p+14) ^-8(p-2)/3(p+14) . 

(A 12 ) 


In the v m 

A'Kjiso 


< v < v c regime, one has 

(vF v (v = 10 18 Hz)) 16 /( p+ 18 ) D 32 ^ P+1S ' ) j“ 16 /h’+ 18 ) n (p-io)/(p+i 8 )£ 3 (p+ 2 )/(p+i 8 ) 

X(1 + ^)-16/(p+18) e -16/(p+18) e -12/(p+18) i/ 8(p-3)/(p+18) ) ( Al3 ) 
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= 


and 


e (1159.5+64.4p)/(p+14) i/ -2(p+18)/3(p+14)^^^ = 1Q 18 Hz ^-16/3(p+14) D -32/3(p+14) f 16/3(p+14) 
Xn -(3p+26)/(3p+42) i —4(p+6)/3(p+14)/ 1 + y )-4(p+18)/3(p+14) ^ + m-(p+2)/3(p+14) 

(A14) 


xe 


16/3(p+14) — 8(p—3)/3(p+14) 


18 


fp 3 = 5.53 x 10 _15 e°* 767p F-A 


(A15) 


A.2. The wind model 


The following derivations follow iGao et al.l ( 20131) ; IChevalier fe Lil (2000 (1 ; iLii fc Zhang 


(2014). 

For the p > 2 case, the forward shock emission can be characterized as 


v m = 5.2 x 10 11 Hz 


P 


p — 1 


(l + z) 1/2 e 1/2 e 2 E 1/2 t~ 3/2 

G ^ *) fc S,-2 fc e,-l- C/ K,iso,52 t d i 


Vc = 1,7 x 10 18 Hz (1 + z)- s/2 (l + H)- 2 el 2 2 BK'!L, 52 4,- 2 _ 1 iy 2 , 


m2 A/2 


F, 


= 1.6 mJy (1 + z) 3/2 D. 


-2 1/2 pl /2 


7? / 4 + 

- C/ K.isn.'i2- rt *,-l t ', 


- 1/2 


28 c B,-2-^K,iso,52 


(A16) 


where A* is the density parameter of the stellar wind medium. 

For p > 2 and in the u > rnax(z/ m , z/ c ) regime, one has 

vF„(v = 10 18 Hz) = F U}mSLX vl /2 v£- 1)/2 v ( x~ p)/2 

= 2.6 x 10" 13 erg s" 1 cm- 2 D^( 1 + z) (p+2)/4 (l + F) _1 /p 

(A17) 


p (p+2)/4 (2—3p)/4 (2-p)/2 
A ^K,iso,52 G /y 18 > 


f P 4 = 6.73 


p — 2 
p — 1 


p-i 


(5.2 x iq- 7 )(p- 2 - 3 )/ 2 , 


(A18) 


Ft< : 


K,iso,52 


z/F iy (/y = 10 18 Hz) 

2.6 x 10~ 13 erg s _1 cm -2 


4/(p+2) 


v 4(l-p)/(p+2) (3p-2)/(p+2) 2(p-2)/(p+2) 
Afc e,-1 L d ms 


b 8/( P + 2) (i + + yj4/(p+2)y-4/(p+2) £ g-r)/(,^2) 

(A19) 


CB.-2 ~ 


//Flz/ = 10 18 Hz) 


i/(p+i) 


1.7 x 10 18 Hz\ (p+2)/(2p+2) 

v c ) ^2.6 x 10 -13 erg s _1 cm -2 

v /i I V'/ -1 f“ 1 /(P+ 1 )A 1 -P)/(P+ 1 ) 4 -(p+2)/(P+ 1 )+p/(p+ 1 ) i .(p-2)/2(p+1) 
X G + J J / p4 e e,-l ^*,-1 l d U 18 


^>28 (P+1) ( 1 + ^)' (P+2)/(P+1) 
(A20) 
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For p > 2 and in the v m < v < v c regime, one has 
vF v (v = 10 18 Hz) = 

= 2,0 X 10 - 13 erg sT 1 + z)^ 1 

(A21) 


v A .(l-3p)/4 (3-p)/2 
I'd rn 8 j 


y K,iso,52 


4/ (p+1) 


xe 


1 /^( 1 / = 10 18 Hz) 

2.0 x 10 -13 erg s _1 cm -2 

4(l-p)/(P+4) A -4/(p+1) , (3p-l)/(p+1) 2(p—3)/(p+1) 


£) 3 /(P +1 )(]_ _|_ z ^J-(P+5)/(p+l) y-4/(p+ 1 ) e -! 


e,-l 


’A ,y ) t 


v 


18 


(A 22 ) 


^1.7 x 10 18 Hzy /2 ^ 


^ B ,-2 — 


z/i^i/ = 10 18 Hz) 


l/(p+l) 


.0 x 10 13 erg s 1 cm 2 

v (1 _j_ V+ — 1 /t-(P+ 2 )/(P+ 1 ) f-l/(p+l) c (!-p)/(p+ 1 )j.p/(P+ 1 ),.(p-3)/2(p+1) 
Xd + l j A*-l /p4 e e,-l hi A 8 


D 2 8 /(p+1) (l + ^)-(p+2)/(p+ 1 ) 
(A23) 


For p < 2 and in the z/ > max(z/ m , y) regime, one has 

E K ,iso = {vF v {v = 10 18 Hz )) 8 /(P+6) jD 1 6/(p+6)y.- 8/( P+6) A (p- 2 )/(p+6) t 

X (1 + Y') 8/(p+6) x (1 + ^)- 6 /(p+ 6 ) e - 8 /(p+ 6 ) I/ 4 (p- 2 )/(p+ 6 )^ (A24) 


f p5 = 405854 x 3 3 P/ 8 5 _ P / 4 e _ 7 p 7 r (4_3p)/8 


2 — p 
p — 1 


(A25) 


£b — 


6.87 x 10- n v- 2/3 (vF u (v = 10 18 Hz)) 8 / 3 (p+ 6 ) D 16 / 3 (p+ 6 ) /p 5 8 / 3 (p+ 6 ) A- ( 3 p+ 26 )/( 3 p+ 18 ) t 2/3 

x (l _|_ y)-4(p+4)/3(p+6)^ _|_ 2; ^-(p+8)/(p+6) e -8/3(p+6) I/ 4(p-2)/3(p+6)^ (A26) 


For p < 2 and in the v m < v < v c regime, one has 

F/ K ,iso = (vF„(v = 10 18 Hz)) 8/(p+4) D 16/(p+4 V p _ 6 8/(p+4) A (p - 10)/(p+4) t (p+8)/(p+4) 

X (1 + ^)- 12 /(p+ 4 ) e -8/(P+ 4 ) e - 6 / ( p+4) z/ 4(P- 3 )/(P+ 4 ), (A27) 

fp6 = 1.72 x 10 22 x 2- 17 P/ 2 3 3p/8 5“ 37p/4 e 13 - 38p 7T (4 “ 3p)/8 (—j) > ( A28 ) 

= (5.69 x lO _16 ) 2 ^ +4 ^ 3 ^ +6 V~ 2< ' p+4 ^ 3 ^ +6 ' ) (z/F iy (z/ = 10 18 Hz)) 8 / 3 (p+ 6 )y) 16 / 3 (p+ 6 )j” 8 / 3 ( p+6 ) 

Xx4 _< ' 3p+26 )/( 3p+18 )t 2 / 3 (l + y)- 4 (P+ 4 )/ 3 (P+6)(l _)_ 2; ^-(p+8)/(p+6) e -8/3(p+6) z/ 4(p-3)/3(p+6) 

(A29) 
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Table 1. The temporal and spectral parameters of 85 GRBs. 


GRB /3q fix a O,2 u Function a.x ,i a a X,2 u Function Aax^* 5 ^fix,0 C 


Grade I 


050408 

0.28 ± 0.33 

1.14 ± 0.14 

0.49 ± 0.01 

1.29 ± 0.11 

3 

BPL 

0.73 ± 0.15 

1.17 ± 0.22 

3 

BPL 

-0.12 ± 0.33 

0.86 ± 0.47 

40.7 

050801 

0.69 ± 0.34 

0.92 ± 0.17 

0.07 ± 0.01 

1.20 ± 0.01 

3 

BPL 

0.24 ± 0.11 

1.18 ± 0.03 

3 

BPL 

-0.02 ± 0.04 

0.23 ± 0.51 

0.2 

050820A 

0.72 ± 0.03 

0.89 ± 0.05 

0.91 ± 0.02 

1.67 ± 0.09 

3 

BPL 

1.12 ± 0.08 

1.89 ± 0.11 

3 

BPL 

0.22 ± 0.20 

0.17 ± 0.08 

2379.0 

050922C 

0.51 ± 0.05 

1.06 ± 0.11 

0.82 ± 0.11 

1.53 ± 0.09 

3 

BPL 

1.04 ± 0.12 

1.71 ± 0.19 

3 

BPL 

0.18 ± 0.28 

0.55 ± 0.16 

8.0 

051028 

0.60 ± 0.00 

0.95 ± 0.15 

0.99 ± 0.06 



SPL 

1.16 ± 0.08 



SPL 

0.17 ± 0.14 

0.35 ± 0.15 


051109A 

0.70 ± 0.05 

0.98 ± 0.08 

0.64 ± 0.08 

1.07 ± 0.12 

3 

BPL 

0.24 ± 0.04 

1.22 ± 0.11 

3 

BPL 

0.15 ± 0.23 

0.28 ± 0.13 

3.5 

060111B 

0.70 ± 0.10 

0.95 ± 0.18 

0.80 ± 0.07 

1.55 ± 0.08 

3 

BPL 

0.90 ± 0.15 

1.59 ± 0.12 

3 

BPL 

0.04 ± 0.20 

0.25 ± 0.28 

7.2 

060206 

0.73 ± 0.05 

1.20 ± 0.31 

0.42 ± 0.09 

1.43 ± 0.10 

3 

BPL 

0.40 ± 0.09 

1.50 ± 0.06 

3 

BPL 

0.07 ± 0.16 

0.47 ± 0.36 

12.5 

060418 

0.78 ± 0.09 

1.10 ± 0.10 

1.23 ± 0.07 


3 

BPL 

1.33 ± 0.06 


3 

BPL 

0.10 ± 0.13 

0.32 ± 0.19 


060512 

0.68 ± 0.05 

1.04 ± 0.10 

0.81 ± 0.05 



SPL 

1.20 ± 0.07 



SPL 

0.39 ± 0.12 

0.36 ± 0.15 


060714 

0.44 ± 0.04 

1.10 ± 0.19 

0.15 ± 0.06 

1.04 ± 0.15 

3 

BPL 

0.48 ± 0.09 

1.34 ± 0.11 

3 

BPL 

0.30 ± 0.26 

0.66 ± 0.23 

5.9 

060729 

0.78 ± 0.03 

1.02 ± 0.04 

0.10 ± 0.05 

1.40 ± 0.15 

3 

BPL 

0.05 ± 0.01 

1.45 ± 0.11 

3 

BPL 

0.05 ± 0.26 

0.24 ± 0.07 

53.0 

060904B 

1.11 ± 0.10 

1.19 ± 0.15 

1.10 ± 0.05 


3 

BPL 

1.41 ± 0.18 


3 

BPL 

0.31 ± 0.23 

0.08 ± 0.25 

2.4 

060912A 

0.60 ± 0.15 

0.62 ± 0.20 

0.94 ± 0.03 



SPL 

1.07 ± 0.02 



SPL 

0.13 ± 0.05 

0.02 ± 0.35 


060927 

0.61 ± 0.05 

0.77 ± 0.20 

1.30 ± 0.10 


3 

BPL 

1.30 ± 0.07 


3 

BPL 

0.00 ± 0.17 

0.16 ± 0.25 

0.9 

061007 

1.02 ± 0.05 

1.00 ± 0.10 

1.62 ± 0.08 


3 

BPL 

1.66 ± 0.07 



SPL 

0.04 ± 0.15 

-0.02 ± 0.15 


061126 

0.82 ± 0.09 

0.85 ± 0.17 

1.29 ± 0.04 


3 

BPL 

1.34 ± 0.05 



SPL 

0.05 ± 0.09 

0.03 ± 0.26 

6.0 ° 

070318 

0.78 ± 0.10 

0.97 ± 0.11 

1.02 ± 0.10 


3 

BPL 

1.03 ± 0.02 



SPL 

0.01 ± 0.12 

0.19 ± 0.21 


070411 

0.75 

1.24 ± 0.22 

0.50 ± 0.08 

1.50 ± 0.11 

3 

BPL 

1.10 ± 0.06 

1.40 ± 0.09 

3 

BPL 

-0.10 ± 0.20 

0.49 ± 0.22 

65.0 

070518 

0.80 

1.20 ± 0.34 

0.70 ± 0.07 

1.80 ± 0.11 

3 

BPL 

0.41 ± 0.06 

1.51 ± 0.09 

3 

BPL 

-0.29 ± 0.20 

0.40 ± 0.34 

40.1 

071025 

0.96 ± 0.14 

1.08 ± 0.11 

1.43 ± 0.06 


3 

BPL 

1.52 ± 0.08 



SPL 

0.09 ± 0.14 

0.12 ± 0.25 


071031 

0.64 ± 0.05 

0.71 ± 0.14 

0.79 ± 0.05 


3 

BPL 

0.82 ± 0.05 



SPL 

0.03 ± 0.10 

0.07 ± 0.19 


080319C 

0.98 ± 0.42 

0.61 ± 0.10 

1.12 ± 0.13 


3 

BPL 

1.33 ± 0.08 



SPL 

0.21 ± 0.21 

-0.37 ± 0.52 


080413A 

0.52 ± 0.37 

1.15 ± 0.24 

1.54 ± 0.05 


3 

BPL 

1.68 ± 0.09 


3 

BPL 

0.14 ± 0.14 

0.63 ± 0.61 

0.3 

080603A 

0.98 ± 0.04 

1.01 ± 0.10 

0.95 ± 0.03 


3 

BPL 

0.96 ± 0.05 



SPL 

0.01 ± 0.08 

0.03 ± 0.14 


080710 

0.80 ± 0.09 

1.00 ± 0.11 

0.39 ± 0.05 

1.32 ± 0.11 

3 

BPL 

0.34 ± 0.04 

1.57 ± 0.14 

3 

BPL 

0.25 ± 0.25 

0.20 ± 0.20 

6.8 

080804 

0.43 

0.82 ± 0.10 

0.87 ± 0.01 



SPL 

1.11 ± 0.01 



SPL 

0.24 ± 0.02 

0.39 ± 0.10 


080913 

0.79 ± 0.03 

1.01 ± 0.23 

0.98 ± 0.02 



SPL 

1.32 ± 0.15 



SPL 

0.34 ± 0.17 

0.22 ± 0.26 


080928 

1.32 ± 0.22 

1.14 ± 0.10 

2.02 ± 0.12 


3 

BPL 

1.81 ± 0.11 


3 

BPL 

-0.21 ± 0.23 

-0.18 ± 0.32 

7.1 

081008 

0.40 ± 0.23 

0.98 ± 0.11 

0.64 ± 0.06 

1.60 ± 0.09 

3 

BPL 

0.87 ± 0.15 

1.68 ± 0.08 

3 

BPL 

0.08 ± 0.17 

0.58 ± 0.34 

9.5 

081203A 

0.60 

1.04 ± 0.10 

1.15 ± 0.07 

1.87 ± 0.13 

3 

BPL 

1.04 ± 0.09 

1.89 ± 0.11 

3 

BPL 

0.02 ± 0.24 

0.44 ± 0.10 

7.1 

090102 

0.74 ± 0.22 

0.79 ± 0.11 

0.20 ± 0.04 

1.16 ± 0.09 

3 

BPL 

0.31 ± 0.05 

1.41 ± 0.09 

3 

BPL 

0.25 ± 0.18 

0.05 ± 0.33 

1.0 

090323 

0.74 ± 0.15 

0.87 ± 0.22 

1.55 ± 0.05 



SPL 

1.62 ± 0.09 



SPL 

0.07 ± 0.14 

0.13 ± 0.37 


090328 

1.19 ± 0.21 

0.90 ± 0.30 

1.84 ± 0.08 



SPL 

1.67 ± 0.11 



SPL 

-0.17 ± 0.19 

-0.29 ± 0.51 


090426 

0.76 ± 0.14 

1.03 ± 0.15 

0.14 ± 0.09 

1.25 ± 0.04 

3 

BPL 

0.13 ± 0.02 

1.04 ± 0.05 

3 

BPL 

-0.21 ± 0.09 

0.27 ± 0.29 

0.2 

090618 

0.50 ± 0.05 

0.92 ± 0.05 

0.76 ± 0.11 

1.53 ± 0.11 

3 

BPL 

0.93 ± 0.09 

1.74 ± 0.10 

3 

BPL 

0.21 ± 0.21 

0.42 ± 0.10 

45.1 

090926A 

0.72 ± 0.17 

0.98 ± 0.15 

1.34 ± 0.05 



SPL 

1.41 ± 0.03 



SPL 

0.07 ± 0.08 

0.26 ± 0.32 


091127 

0.18 

0.68 ± 0.11 

0.55 ± 0.11 

1.50 ± 0.11 

3 

BPL 

0.96 ± 0.05 

1.59 ± 0.12 

3 

BPL 

0.09 ± 0.23 

0.50 ± 0.11 

35.3 

100418A 

0.98 ± 0.09 

1.04 ± 0.29 

0.11 ± 0.01 

1.60 ± 0.10 

3 

BPL 

-0.12 ± 0.03 

1.57 ± 0.11 

3 

BPL 

-0.03 ± 0.21 

0.06 ± 0.38 

90.1 

100901A 

0.52 ± 0.10 

1.00 ± 0.30 

1.42 ± 0.02 


3 

BPL 

1.41 ± 0.02 


3 

BPL 

-0.01 ± 0.04 

0.48 ± 0.40 

29.8 

101024A 

0.70 ± 0.40 

0.82 ± 0.13 

0.01 ± 0.05 

1.07 ± 0.08 

3 

BPL 

-0.09 ± 0.07 

1.37 ± 0.10 

3 

BPL 

0.30 ± 0.18 

0.12 ± 0.53 

1.0 

120326A 

0.75 ± 0.08 

0.77 ± 0.06 

1.52 ± 0.10 


3 

BPL 

1.69 ± 0.09 


3 

BPL 

0.17 ± 0.19 

0.02 ± 0.14 

35.50 

130427A 

0.69 ± 0.01 

0.68 ± 0.16 

1.02 ± 0.09 

1.84 ± 0.11 

3 

BPL 

1.09 ± 0.07 

1.63 ± 0.09 

3 

BPL 

-0.21 ± 0.20 

-0.01 ± 0.17 

127.5 


Grade II 


051111 

0.78 ± 0.07 

1.24 ± 0.17 

1.56 ± 0.11 

3 

BPL 

1.60 ± 0.12 

SPL 

0.04 ± 0.23 

0.46 ± 0.24 

3.0 ° 

090313 

0.74 ± 0.40 

1.08 ± 0.17 

1.55 ± 0.12 

3 

BPL 

1.67 ± 0.10 

SPL 

0.12 ± 0.22 

0.34 ± 0.57 

20.5 ° 


Grade III 


050319 

0.74 ± 0.42 

1.01 ± 0.07 

0.39 ± 0.06 

1.02 ± 0.04 

3 

BPL 

0.58 ± 0.07 

1.72 ± 0.11 

3 

BPL 

0.70 ± 0.15 

0.27 ± 0.49 

55.0 

050401 

0.50 ± 0.20 

0.79 ± 0.13 

0.50 ± 0.08 

0.89 ± 0.08 

3 

SPL 

0.76 ± 0.05 

1.67 ± 0.11 

3 

BPL 

0.78 ± 0.19 

0.29 ± 0.33 

4.3 X 

050416A 

1.30 

1.07 ± 0.11 

0.26 ± 0.07 

1.12 ± 0.12 

3 

BPL 

0.66 ± 0.12 

0.99 ± 0.02 

3 

BPL 

-0.13 ± 0.14 

-0.23 ± 0.11 

11.0 

050525A 

0.52 ± 0.08 

1.09 ± 0.16 

1.46 ± 0.07 


3 

BPL 

1.57 ± 0.04 


3 

BPL 

0.11 ± 0.11 

0.57 ± 0.24 

4.2 ° 

050603 

0.20 ± 0.10 

1.02 ± 0.13 

1.70 ± 0.13 



SPL 

1.71 ± 0.05 



SPL 

0.01 ± 0.18 

0.82 ± 0.23 


050721 

1.16 ± 0.35 

0.85 ± 0.22 

0.60 ± 0.03 



SPL 

1.01 ± 0.08 



SPL 

0.41 ± 0.11 

-0.31 ± 0.57 


050730 

0.52 ± 0.05 

1.62 ± 0.04 

0.48 ± 0.05 

1.47 ± 0.06 

3 

BPL 

0.45 ± 0.13 

2.64 ± 0.20 

3 

BPL 

1.17 ± 0.26 

1.10 ± 0.09 

90.1 

051221A 

0.64 ± 0.05 

1.06 ± 0.14 

0.34 ± 0.07 

1.24 ± 0.04 

3 

BPL 

0.35 ± 0.08 

1.34 ± 0.04 

3 

BPL 

0.10 ± 0.08 

0.42 ± 0.19 

25.1 

060210 

0.37 ± 0.08 

1.08 ± 0.08 

0.53 ± 0.05 

1.77 ± 0.14 

3 

BPL 

0.53 ± 0.06 

1.30 ± 0.12 

3 

BPL 

-0.47 ± 0.26 

0.71 ± 0.16 

5.0 

060526 

0.51 ± 0.32 

0.90 ± 0.11 

0.56 ± 0.10 

1.93 ± 0.10 

1 

BPL 

0.67 ± 0.08 

2.06 ± 0.15 

3 

BPL 

0.13 ± 0.25 

0.39 ± 0.43 

50.1 

060605 

1.06 

1.02 ± 0.09 

0.13 ± 0.09 

2.64 ± 0.15 

3 

BPL 

0.55 ± 0.08 

2.82 ± 0.15 

3 

BPL 

0.18 ± 0.30 

-0.04 ± 0.09 

15.0 

060614 

0.47 ± 0.04 

0.90 ± 0.09 

-0.35 ± 0.04 

1.90 ± 0.12 

3 

BPL 

0.11 ± 0.10 

1.97 ± 0.15 

3 

BPL 

0.07 ± 0.27 

0.43 ± 0.13 

44.0 

060906 

0.56 ± 0.02 

1.08 ± 0.17 

1.45 ± 0.21 


3 

BPL 

1.25 ± 0.11 


3 

BPL 

-0.20 ± 0.32 

0.52 ± 0.19 

1.3 

060906 

0.56 ± 0.02 

1.08 ± 0.17 

1.44 ± 0.20 


3 

BPL 

1.90 ± 0.11 


3 

BPL 

0.46 ± 0.31 

0.52 ± 0.19 

10.5 


o H o 0 
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Table 1—Continued 


GRB 

/5o 

Px 


Q 

o,i a 

Q 

O 

2 

UJ 

Function 

cxx 

a 

1 

ocx,2 

w 

Function 

A “X,O b 

A/3 X ,O c 

+ d 

060908 

0.24 

± 0.20 

1.13 

± 

0.19 

0.71 

± 

0.04 

1.11 

± 

0.09 

3 

BPL 

0.54 

± 

0.05 

1.66 

± 

0.05 

3 

BPL 

0.55 

± 

0.14 

0.89 

± 

0.39 

1.1 

070110 

0.55 

± 0.04 

1.08 

± 

0.11 

0.16 

± 

0.06 

1.67 

± 

0.12 

3 

BPL 

0.30 

± 

0.10 

5.10 

± 

0.30 

3 

BPL 

3.43 

± 

0.42 

0.53 

± 

0.15 

20.3 

070125 

0.59 

± 0.10 

1.03 

± 

0.20 

2.96 

± 

0.05 





BPL 

2.12 

± 

0.04 




3 

BPL 

-0.84 

± 

0.09 

0.44 

± 

0.30 

101.1 

070306 

0.70 

± 0.10 

0.95 

± 

0.09 

1.22 

± 

0.11 




3 

BPL 

2.03 

± 

0.04 




3 

BPL 

0.81 

± 

0.15 

0.25 

± 

0.19 

36.9 

070311 

1.00 

± 0.20 

1.00 

± 

0.24 

0.73 

± 

0.02 





SPL 

1.09 

± 

0.06 





SPL 




0.00 

± 

0.44 


070419A 

0.48 

± 0.48 

1.20 

± 

0.30 

1.28 

± 

0.04 




3 

BPL 

0.60 

± 

0.02 




3 

BPL 

-0.68 

± 

0.06 

0.72 

± 

0.78 


071010A 

0.61 

± 0.12 

1.30 

± 

0.40 

2.19 

± 

0.08 




3 

BPL 

1.89 

± 

0.07 




3 

BPL 

-0.30 

± 

0.15 

0.69 

± 

0.52 

70.1 

071112C 

0.63 

± 0.29 

0.67 

± 

0.13 

0.30 

± 

0.07 

0.95 

± 

0.11 

3 

BPL 

0.50 

± 

0.08 

1.49 

± 

0.11 

3 

BPL 

0.54 

± 

0.22 

0.04 

± 

0.42 

1.5 

080310 

0.42 

± 0.12 

0.95 

± 

0.18 

0.11 

± 

0.01 

1.24 

± 

0.02 

3 

BPL 

0.03 

± 

0.06 

1.24 

± 

0.08 

3 

BPL 

0.00 

± 

0.10 

0.53 

± 

0.30 

5.1 

080319A 

0.77 

± 0.02 

0.89 

± 

0.10 

0.65 

± 

0.07 




3 

BPL 

0.94 

± 

0.05 






0.29 

± 

0.12 

0.12 

± 

0.12 


080319B 

0.51 

± 0.26 

0.81 

± 

0.07 

0.93 

± 

0.11 

1.60 

± 

0.12 

3 

BPL 

0.73 

± 

0.05 

2.73 

± 

0.16 

3 

BPL 

1.13 

± 

0.28 

0.30 

± 

0.33 

3.0 

080319B 

0.51 

± 0.26 

0.81 

± 

0.07 

1.03 

± 

0.10 

1.69 

± 

0.09 

3 

BPL 

1.43 

± 

0.09 

2.19 

± 

0.11 

3 

BPL 

0.50 

± 

0.20 

0.30 

± 

0.33 

690.7 

080413B 

0.25 

± 0.07 

0.94 

± 

0.07 

0.31 

± 

0.15 

1.89 

± 

0.22 

3 

BPL 

0.92 

± 

0.16 

1.91 

± 

0.23 

3 

BPL 

0.02 

± 

0.45 

0.69 

± 

0.14 

148.5 

080721 

0.68 

± 0.02 

0.94 

± 

0.06 

1.17 

± 

0.03 

1.31 

± 

0.05 

3 

BPL 

0.81 

± 

0.01 

1.65 

± 

0.07 

3 

BPL 

0.34 

± 

0.12 

0.26 

± 

0.08 

3.1 

090510 

0.85 

± 0.05 

0.75 

± 

0.12 

0.84 

± 

0.05 




3 

BPL 

2.27 

± 

0.06 




3 

BPL 

1.43 

± 

0.11 

-0.10 

± 

0.17 

1.5 

090812 

0.36 

0.89 

± 

0.14 

1.27 

± 

0.05 




3 

BPL 

1.22 

± 

0.09 





SPL 

-0.05 

± 

0.14 

0.53 

± 

0.14 


091029 

0.49 ± 0.12 

1.12 

± 

0.08 

0.48 

± 

0.06 

1.34 

± 

0.09 

3 

BPL 

0.32 

± 

0.09 

1.35 

± 

0.08 

3 

BPL 

0.01 

± 

0.17 

0.63 

± 

0.20 

20.8 

100219A 

0.56 

0.69 

± 

0.23 

0.74 

± 

0.08 

1.91 

± 

0.12 

3 

BPL 

0.54 

± 

0.07 

1.65 

± 

0.15 

3 

BPL 

-0.26 

± 

0.27 

0.13 

± 

0.23 

1.8 

100219A 

0.56 

0.69 

± 

0.23 

2.21 

± 

0.13 




3 

BPL 

2.51 

± 

0.16 




3 

BPL 

0.30 

± 

0.29 

0.13 

± 

0.23 

20.5 

110205A 

0.49 ± 0.08 

0.78 

± 

0.06 

1.51 

± 

0.08 




3 

SPL 

1.59 

± 

0.02 





SPL 

0.08 

± 

0.10 

0.29 

± 

0.14 


110918A 

0.42 

± 0.18 

0.89 

± 

0.30 

1.65 

± 

0.07 





SPL 

1.61 

± 

0.12 





SPL 

-0.04 

± 

0.19 

0.47 

± 

0.48 


120729A 

1.00 ± 0.10 

0.80 

± 

0.17 

0.94 

± 

0.05 

2.27 

± 

0.09 

3 

BPL 

1.09 

± 

0.06 

2.40 

± 

0.16 

3 

BPL 

0.13 

± 

0.25 

-0.20 

± 

0.27 

6.61 

120815A 

0.78 ± 0.01 

0.72 

± 

0.11 

0.63 

± 

0.04 




3 

BPL 

0.86 

± 

0.06 





SPL 

0.23 

± 

0.10 

-0.06 

± 

0.12 



Grade IV 


070611 

0.73 ± 0.00 

0.83 ± 0.31 

0.58 ± 0.04 

3 

BPL 

1.34 ± 0.12 

3 

BPL 

0.76 ± 0.16 

0.10 ± 0.31 

071003 

0.35 ± 0.11 

0.91 ± 0.12 

1.62 ± 0.11 

3 

BPL 

1.63 ± 0.02 


SPL 

0.01 ± 0.13 

0.56 ± 0.23 

120711A 

0.52 ± 0.02 

0.81 ± 0.10 

0.96 ± 0.04 

3 

BPL 

1.64 ± 0.05 


SPL 

0.68 ± 0.09 

0.29 ± 0.12 


Grade V 


060607A 

070208 

070420 


0.72 ± 0.27 0.62 ± 0.06 -0.93 

0.68 1.20 ± 0.20 0.49 ± 0.06 

-1.43 ± 0.04 


4.60 

0.90 ± 0.08 


3 

3 

3 


BPL 

BPL 

BPL 


0.36 ± 0.03 
0.43 ± 0.03 
0.12 ± 0.01 


3.10 ± 0.12 
1.78 ± 0.13 
1.46 ± 0.05 


3 

3 

3 


BPL 

BPL 

BPL 


-0.10 ± 0.33 9.5 x 
0.52 ± 0.20 9.0 X 

3.0 X 


a For single power-law (SPL) decay lightcurves (as described in section 3), the decay indices are also denoted as cki; 
bA “X.O = “X,2 - “ 0 , 2 ; 

°A/3 x ,o = Px - Po; 

[11 units of ks. The symbols “X” and “O” denote the X-ray and optical bands, respectively. 

















Table 2. The temporal decay index a and spectral index f3 in different afterglow models. 


CMB 

Spectral regime 

Pip) 

a{p)/a{p, q) 

q) 

a(p)/a(p,q) 

«(/3)/a(/3, <?) 





p> 2 


1 < p < 2 

Adiabatic deceleration without energy injection 

ISM 

Wind 

Vm < V < 

V > V c 

Urn <V<V C 

V > 1S C 

P-1 

2 

V 

A 

2 

V 

_2_ 

3p 4 -2 

3p-l 

3 p-2 

_4_ 

a= M 

LX. 2 

rv 3 ' 0_1 

a = A 

a - ^ 

___2_ 

3pH?10 

P+8 

pA 

_8_ 

a= M±2 

«=# 
a= ^ 

a= M±6 

___8_ 


Adiabatic deceleration with energy injection 


ISM 

Wind 

Vm < V < 

V > V c 

I'm < V < U c 

V > I'c 

P-1 

(2p-6) + (p+3jq 

A (a 11 1 (2+9) ^ 

lJ-18,-p( g +2) 

199-10 , (2+91/3 

2 

V 

2 

p-1 

2 

V 

_2_ 

(2p-4)+(p+2)q 

(2p-2) + (p+l)q 

(2p-4)+(p+2)q 
_4_ 

Q w A / i o 

_ <7-2 , (2+q)(3 

a=i + ^ 

a _ ,-2 + (2+,W> 

16 x 

14q+p(q+2)—4 

,16 x 
4+(p+4)q 

(6 +p)q 

__8_ 

01 16 + 8 
.. _ 79-2 , (2+9)/3 
^ 8 ' 8 
_ 69±4 §q 

8 T 4 
„ - <0+3)9 

___4_ 

Post jet break phase without energy injection 

ISM 

Wind 

I'm < V < I'c 

V > I'c 

I'm < V < I'c 

I' > I'c 

p-1 

2 

P 

A 

2 

V 

_2_ 

eo| col 

_ -4 _ 

307+6] 

16 

3p+22 

P+ 6 12 

pAo 

_8_ 

_ 3(2/l+7) 

16 

a = 3 < 3 + 11 
„ - 2/3+13 

a^-,4 5 

Post jet break phase with energy injection 

ISM 


p-1 

p(q+2)—4(1 —q) 

5q —2 . (2+q)P 

22q — 4+p(q+2) 

119-2 , (2+91/1 








v > v c 

V 

2 

3q — 2+p(q+2) 

4 

_ 3q—2+2/3(q+2) 

' ' 4 

1 8( j+4+p(9+2) 

16 

_ 9 9 +2+/3(9+2) 

8 

Wind 

I'm < V < I'c 

p-1 

3q — 2+p(q+2) 

i (2+q)/3 
a = q + v 

pq+Sq+4 

„ - 1 | (2/3+9)9 


V > I'c 

P 

p(9+2)-4(l- 9 ) 

_ /3(q+2)-2(l- g ) 

(p+10)9 

a 1 (/3+5)9 
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Table 3. 


The A a?x,o 


and A /3 x ,o values in different afterglow models, p > 2. 


p > 2 

Same regime a 

Different 

regimes** 

Grey 

zone c 

ISM,wind 

ISM 

wind 

ISM 

wind 


o o 

II II 

0 0 

* x 

3 

< < 

A /3x,o=i 

A “X.O=3 

A /3 x,o = * 
A “X,0=-3 

A /3x,O = (0,*) 
A “X,O=(0,|) 

A /3x,O = (0.i) 
A “X.O = (-3.°) 



Energy i 

njection break 


A ^x,o = 0 

A “l,X,O =0 

Aa 2,X,O= 0 

A /3 X,0=3 

A “i,x,o = ^ir- 

A “2,X,0 = 3 

a /5x.O = 5 

A “l,X.O = ~4 +<i 
A “2,X,0=-3 

A /3x,O = (0,i) 
A “l,X.O = 

A “2,X,O = (0,z) 

A /3x,O = (0.4) 

At *l,X,0 = 

Aa 2,X,0 = (-3-°) 


A/3x,O=0 

A “l,X,O =0 

Aa 2,X,O =0 

a /3x.O=5 
A “1,X,0 = 3 
A “2,X.O = T 

a /3x,o = 5 

A “l,X,0=-3 

Aa 2,X,0=-3 

A /3x,0 = ( 0,^) 

A “l,X,O = (0.f ) 
A “2,X,O = (0,t) 

A /?x,o=( o.4) 
Aa i,x,o=(-f -0) 
Aa 2,X,0 = (-|.°) 

Energy injection with jet break 

A /3x.O = 0 
A “i,x,o=° 

Aa 2,X,O=0 

A /3 X.O=i 

A “i,x,o = 
a “2,x,o = 

a /5x.O = 5 

A “l,X.O = ~4 +<i 
A “2,X,0 = A +g 

A /3x,O = (0 

A “l,X.O = ( 0 
a “2,X.O = (0. T 1 ) 

A /3x,O = (0.^) 
A “i,x,o = (0,\^) 
a “2,x,o = (0, 2-r 1 ) 


a Same regime: X-ray and optical bands are in the same spectral regime, regime I (i/ > u c ) or regime II 
(yrn < v < u c ). In this table, A/?x O = fix ~ fiOi and Accx O = a X — ^O- The subsripts “1” and “2” denote 
the pre- and post-break segments for the BPL lightcurves, respectively. 

^Different regimes: X-ray band in regime I (u > i/ c ), and optical band in regime II (i/ m < v < Vc)\ 
c Grey zone: One band or two bands in the grey zone regime I— II. 
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Tablc 4. 


The A ax,o and A/3 Y ,o values in different afterglow models, 1 < p < 2. 


1 < V < 2 

Same regime a 

Different 

regimes** 

Grey 

zone c 

ISM,wind 

ISM 

wind 

ISM 

wind 


o o 

II II 

0 0 

* x 

<0. 3 

<1 <1 

A /3x,o=i 

A “X.O=! 

A /3x,o = * 
Aa x,o = i 

A /3x,O = (0,*) 
A “X,O=(0,|) 

A /3x,O = (0,4) 
A “X,0 = (ol) 



Energy i 

njection break 


A ^X ,0 = 0 

A /3x,0 = + 
A “i,x,o = 

a /3x,0 = 5 

A /3x,O = (0,i) 

A /3x,O = (0.4) 
A «1,X,0 = 

A “1 ,x,o=° 

a “i,x,o = A +g 

A “1,X,0 = (0. Aj- 1 ) 

Aci 2 ,X ,0 = 0 

^ Q 2 ,X, 0 = I 

^ a 2 ,X ,0 = I 

a “2,X,O=(0.1) 

Aa 2,X,O = (0,z) 


A/3x,O=0 

A / 3 x.O =5 

A /5x.O = 5 

A /?x,o=( o,f) 

A /?x,o=( o, 2 ) 

A “l,X,O =0 

Aa 2 ,X,O =0 

A “l,x ,0 = 5 

Aa 2,X,0 = l 

Ac *1,X,0 = t 
A « 2 ,X, 0 =-I 

A “l,X,O = ( 0 d) 
A “2,X,O = (0,|) 

Aa i,x,o=(°»i) 

Aq: 2,X,O = (--Z»0) 

Energy injection with jet break 

A /3x,O=0 

a /3x,0=3 

A /3x.O = 5 

A /3x,O = (0,|) 

A /3x,O = (0 2 ) 
Ac *i,x,o = (A^a.o) 

A “l,X,O=0 

A “1,X,0 = Ay 3 ' 

A “l,X.O = — 4 +g 

A “i,x,o = 

Ao:2,X,0=° 

A “2,X,0 = T 1 

A “2,X,0 = ~ 4 +g 

A “2,X,0 = (0. T 2 ) 

a «2,x,o = (A^ta.o) 


a Same regime: X-ray and optical bands are in the same spectral regime, regime I (i/ > u c ) or regime II 
(yrn < v < u c ). In this table, A/?x O = fix ~ fiOi and Accx O = a X — ^O- The subsripts “1” and “2” denote 
the pre- and post-break segments for the BPL lightcurves, respectively. 

^Different regimes: X-ray band in regime I (u > i/ c ), and optical band in regime II (i/ m < u < v c )\ 
c Grey zone: One band or two bands in the grey zone regime I— II. 
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Table 5. The criteria for GRB grades 


Grades 

Light curve behavior 

Closure-relation 

Fraction 

Grade I 

Achromatic break or SPL decay in both bands 

Yes 

43/85 

Grade II 

Break missing in one band, consistent with being achromatic 

Yes 

2/85 

Grade III 

Achromatic break or SPL decay in both bands 

No 

34/85 

Grade IV 

Break missing in one band, consistent with being achromatic 

No 

3/85 

Grade VI 

Chromatic 

No 

3/85 






Table 6. The Gold sample GRBs and their derived parameters 


GRB 

z 

Optical a 

Xray a 

P 


q 


b 

T7 1 c 

■ C/ 'y ,i so 

■®K,end 

B K ,in C 

^K,dec 

V-y ,dec 

"H-f ,end 

type e 


050408 

1.24 

windll 

windl-windll 

2.12 ± 0.12 

0.25 

± 

0.17 


11.76 

3410.39 ± 76.12 

3386.60 ± 75.06 

23.80 ± 6.84 

33 ± 10 

0.3 ± 0.0 

1 

050801 

1.56 

ISMII 

ISMII 

2.78 ± 0.28 

0.22 

± 

0.15 

2.7E-05 

0.52 ± 0.11 

124.89 ± 17.39 

81.83 ± 20.10 

43.06 ± 17.65 

1 ± 1 

0.4 ± 0.1 

1 

051109A 

2.35 

windl-windll 

windl-windll 

2.4 ± 0.1 

0.22 

± 

0.06 


9.15 ± 6.83 

14675.86 ± 763.28 

14060.42 ± 730.61 

615.44 ± 38.42 

1 ± 1 

0.1 ± 0.1 

1 

060206 

4.05 

windl-windll 

windl 

2.46 ± 0.1 

0.41 

± 

0.06 


4.78 ± 20.79 

386.76 ± 93.02 

344.28 ± 83.18 

42.47 ± 15.11 

10 ± 44 

1 ± 5 

1 

060714 

2.71 

windl-windll 

windl-windll 

1.89 ± 0.07 

0.19 

± 

0.05 


18.22 ± 2.53 

250.46 ± 248.11 

239.83 ± 237.24 

10.62 ± 11.25 

63 ± 67 

7 ± 7 

1 

060729 

0.54 

windl-windll 

windl-windll 

2.04 ± 0.08 

0.06 

± 

0.01 


0.75 ± 0.07 

2312.14 ± 33.54 

2304.60 ± 33.43 

7.54 ± 0.12 

9±1 

0.03 ± 0.01 

1 

070411 

2.95 

windl-windll 

windl-windll 

2.48 ± 0.02 

0.86 

± 

0.09 


17.06 ± 1.84 

24456.57 ± 256.42 

14460.59 ± 794.78 

9995.98 ± 782.91 

0.2 ± 0.1 

0.1 ± 0.1 

1 

070518 

1.16 

windl-windll 

windl-windll 

2.4 ± 0.2 

0.24 

± 

0.08 


0.27 ± 0.13 

917.47 ± 22.99 

911.31 ± 22.81 

6.16 ± 0.52 

4 ± 2 

0.03 ± 0.01 

1 

080710 

0.85 

ISMII 

ISMII 

2.78 

0.32 

9.6E-06 

1.34 ± 0.32 

424.26 ± 40.68 

393.99 ± 37.78 

30.27 ± 2.90 

4 ± 1 

0.3 ± 0.1 

1 

090102 

1.55 

windl-windll 

windl-windll 

1.58 ± 0.22 

0.33 

± 

0.05 


22.74 ± 2.12 

1750.25 ± 1149.59 

1490.68 ± 1009.33 

259.57 ± 141.45 

8 ± 4 

1 ± 1 

1 

090426 

2.61 

windl-windll 

windl 

2.06 ± 0.3 

0.13 

± 

0.02 


0.83 ± 0.28 

2.77 ± 3.33 

1.80 ± 2.17 

0.97 ± 1.17 

46 ± 58 

23 ± 29 

1 

100418A 

0.62 

ISMII 

ISMII 

2.96 ± 0.18 

0.05 

± 

0.02 

2.9E-06 

0.14 ± 0.02 

1754.24 ± 403.22 

1754.06 ± 403.18 

0.18 ± 0.05 

43 ± 14 

0.01 ± 0.01 

1 

101024A 


ISMII 

ISMII 

2.64 ± 0.26 

0.06 

± 

0.11 








1 

050820A 

2.61 

ISMI-ISMII 

ISMI-ISMII 

1.78 ± 0.1 





114.67 ± 33.59 



4117.25 ± 2462.50 

3 ± 2 


2 

050922C 

2.20 

ISMII 

ISMI 

2.06 ± 0.05 





9.93 ± 1.06 



118.25 ± 17.29 

8 ± 1 


2 

060111B 


ISMI 

ISMI 

1.57 ± 0.03 





11.00 ± 5.00 






2 

081008 

1.97 

ISMII 

ISMI 

1.96 ± 0.18 





9.30 ± 1.91 



253.98 ± 93.28 

4 ± 3 


2 

081203A 

2.10 

ISMI-ISMII 

ISMI-ISMII 

2.2 ± 





34.71 ± 17.12 



594.07 ± 14.99 

6 ± 3 


2 

090618 

0.54 

ISMII 

ISMI 

1.92 ± 0.02 





25.30 ± 0.00 



89.10 ± 14.81 

22 ± 4 


2 

091127 

0.49 

ISMII 

ISMI 

1.36 ± 





1.52 ± 0.08 



263.04 ± 22.94 

1 ± 0.1 


2 

130427A 

0.34 

ISMII 

ISMII 

2.38 ± 0.02 




9.1E-06 

81.00 



1665.96 ± 62.14 

5 ± 0.2 


2 

051028 

3.70 

ISMI-SIMII 

ISMI-ISMII 

1.9 ± 0.1 





11.94 ± 1.87 



2515.02 ± 1557.52 

0.5 ± 0.3 


4 

060418 

1.49 

ISMII 

ISMI-ISMII 

2.58 ± 0.18 





12.45 ± 4.44 



386.85 ± 42.95 

3 ± 1 


4 

060512 

0.44 

ISMI-SIMII 

ISMI-ISMII 

1.94 ± 0.06 





0.02 ± 0.01 



5.70 ± 0.12 

0.3 ± 0.2 


4 

060904B 

0.70 

ISMI-SIMII 

ISMI-ISMII 

3.25 ± 0.17 





0.79 ± 0.56 



750.73 ± 388.34 

0.1 ± 0.1 


4 

060912A 

0.94 

ISMII 

ISMII 

2.2 ± 0.3 




8.7E-05 

1.12 ± 0.08 



46.23 ± 1.16 

2 ± 0.2 


4 

060927 

5.46 

windll 

windl-windll 

2.22 ± 0.1 





10.37 ± 3.03 



37039.49 ± 2373.94 

0.03 ± 0.01 


4 

061007 

1.26 

ISMII 

ISMII 

2.04 ± 0.1 




7.3E-05 

99.81 ± 7.22 



1085.73 ± 18.15 

8 ± 1 


4 

061126 

1.16 

ISMII 

ISMII 

2.64 ± 0.18 




1.0E-05 

13.26 ± 1.90 



1282.64 ± 177.08 

1 ± 0.1 


4 

070318 

0.84 

ISMII-ISMII 

ISMI-ISMII 

2.74 ± 0.02 





1.39 ± 0.35 



277.28 ± 22.56 

0.5 ± 0.1 


4 

071025 

1.55 

ISMII 

ISMII 

3.07 ± 0.13 




9.2E-06 

85.42 ± 12.42 



1394.60 ± 851.47 

6±4 


4 

071031 

2.69 

ISMI 

ISMI 

2.01 ± 0.27 





4.73 ± 9.26 



68.53 ± 38.27 

6 ± 5 


4 

080319C 

1.95 

windl-II 

windl-windll 

2.22 ± 0.2 





23.54 ± 1.47 



12314.17 ± 3183.23 

0.2 ± 0.1 


4 

080413A 

2.43 

windll 

windl-windll 

2.3 ± 0.48 





14.09 ± 7.30 



249.98 ± 71.42 

5 ± 3 


4 

080603A 

1.68 

ISMI 

ISMI 

2.96 ± 0.08 





2.20 ± 0.80 



3577.17 ± 272.70 

0.1 ± 0.1 


4 

080804 

2.20 

windl-windll 

windl-windll 

1.86 





24.61 ± 4.79 



143.09 ± 16.57 

15 ± 3 


4 

080913 

6.70 

ISMI-SIMII 

ISMI-ISMII 

2.79 ± 0.27 





8.44 ± 1.55 



354.60 ± 174.38 

2 ± 1 


4 

080928 

1.69 

ISMII 

ISMII 

2.44 ± 0.04 




5.7E-05 

6.30 ± 0.75 



548.00 ± 26.66 

1 ± 0.1 


4 

090323 

3.57 

windll 

windll 

2.65 ± 0.13 




9.5E-04 

372.38 ± 16.86 



209526.73 ± 18842.90 

0.2 ± 0.1 


4 

090328 

0.74 

ISMII 

ISMII 

3.19 ± 0.21 




3.4E-06 

19.03 



3861.76 ± 2039.23 

0.5 ± 0.3 


4 

090926A 

2.11 

ISMII 

ISMII 

2.29 ± 0.19 





185.13 ± 9.10 



1859.10 ± 111.85 

9 ± 1 


4 

100901A 

1.41 

windll 

windl-windll 

2.1 ± 0.1 





2.95 ± 0.63 



27847.74 ± 871.45 

0.01 ± 0.01 


4 

120326A 

1.80 

windll 

windll 

2.51 ± 0.09 




5.2E-04 

3.18 ± 0.40 



41757.87 ± 1069.46 

0.01 ± 0.01 


4 

Grade II 

051111 

1.55 

windll 

windl-windll 

2.52 ± 0.14 





10.79 ± 3.07 



3793.13 ± 273.16 

0.3 ± 0.1 


4 

090313 

3.38 

windll 

windl-windll 

2.54 ± 0.72 




5.0E-06 

13.02 ± 2.94 



173028.58 ± 189747.00 

0.01 ± 0.01 


4 

















a ISMI: the ISM model in the spectral regime I (u > u c )\ ISMII: the ISM model in the spectral regime II (y-m < v < I'c); windl: the wind model in the spectral regime I; windll: the wind model 
in the spectral regime II: ISMI-ISMII: the ISM model in the grey zone (between u > u c and i/ m < v < Vc)', windl-windll: the wind model in the grey zone. 

^Upper limit, calculate with the X-ray data in the spectral regime II (Um < v < I'c); 

c In units of 10^ erg. £?K,end is the kinetic energy at the end of energy injection (tend); ^K,in is the injected kinetic energy during the energy injection phase; and £?K,dec is the kinetic energy 
at the fireball deceleration time (tdec)> 

^In units of %; 77^dec i s the radiative efficiency calculated using £^K,decj and Tj^dec is the radiative efficiency calculated using ^K,end 
e l: energy injection break; 2: jet break; 3: jet break with energy injection; 4: SPL decay. 
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Table 7. Parameters of the jet break sample. 


GRB 

d °j 

log(E 7 /erg) 

l°g(E K /erg) 

050820A 

4.5 ± 3.0 

51.54 ±°;“ 

53.05 1°; 2 9 

050922C 

060111B a 

1.8 ± 0.3 

49.69 1°;? o 

50.13 mi 

081008 

1.3 ± 0.4 

49.4i mt 

50.54 1q29 

081203A 

1.0 ± 0.6 

49-76 1°; 2 6 

50.72 mi 

090618 

3.5 ± 1.3 

50.68 +_° Q l\ 

50.84 mi 

091127 

2.7 ± 0.4 

49 22 + 0 07 
t * t ’.** _q 08 

50 73 + 0 07 
ou ‘° -0.09 

130427A 

3.8 ± 0.3 

51 25 +0 - 04 
OI.ZO _o 04 

51 54 +0-04 

01.04 _ 0- 04 


no redshift z available. 
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Time since trigger (s) 



Time since trigger (s) 




Fig. 1.— The X-ray and optical light curves, as well as the fitting results (blue dotted-dashed 
lines) for the GRBs in the Grade I sample. If an achromatic break exists, the achromatic 
break time is shown by a purple vertical dashed line. 
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Fig. [T]—Continued 
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Fig. [T]—Continued 
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Fig. [T]—Continued 
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Fig. [T]—Continued 
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Fig. 2 


Same as Figured] but for the Grade II sample 
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Fig. 3.— Same as Figure [U but for the Grade III sample 






















- 52 - 



Time since trigger (s) 
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Fig. [3]—Continued 
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Fig. [3]—Continued 
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Fig. [3]— Continued 
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10' 1 
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Fig. 4 


Same as Figured! but for the Grade IV sample 
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Fig. 5.— Same as Figured! but for the Grade V sample 
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Fig. 6.— Typical light curve behaviors: (a) The canonical X-ray light curve, reproduced 
from Zhang et al (2006), with four characteristic temporal segments marked; (b) The energy 
injection break case, with a transition from the shallow decay phase (segment II) to the 
normal decay phase (segment III); (c) The jet break cases without (solid) or with (dashed) 
energy injection. For the former, it is a transition from the normal decay phase (segment 
III) to the post-jet-break phase (segment IV); for the latter, both segments have a shallower 
decay slope; (d) the single power-law (SPL) case. The steep decay phase in the X-ray light 
curve and the early rising phase in the optical light curve are not included in the analysis. 












Fig. 7.— The confrontation of the data with model predictions in the A/3x,o — Aax,o plane. 
The (A f3x,o, A ax,o) model predictions for the ISM and wind medium in different spectral 
regimes are denoted as large grey symbols: squares for ISM and circles for wind, respectively. 
The red dashed lines define the (A (3x,ch Aa^ 0 ) range for the grey zones. The SPL, energy 
injection, and jet break samples are presented by square, circle, and star, respectively. The 
ISM and wind samples are marked with filled and open symbols, respectively. The X-ray 
and optical bands in the same spectral regime, different spectral regimes, and grey zone are 
marked with black, red, and blue, respectively, (a): the entire Gold sample GRBs; (b)and 
(c): the ISM models with p > 2 and 1 < p < 2, respectively; (d) and (e): ththe wind models 
with p > 2 and 1 < p < 2, respectively; (f): the silver sample GRBs (denoted in triangles). 
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Fig. 8.— The distributions of the temporal decay index a for various sub-samples in the 
Gold sample and their best Gaussian fits (dashed curves): (a) all the BPL GRBs in the 
Gold sample, with best fits ao,i = 0-49 ± 0.45, «o ,2 = 1-44 ± 0.39, ax,i = 0.58 ± 0.63, and 
«x ,2 = 1.50±0.27; (b) all the SPL GRBs in the Gold sample, with best fits «o = 1.26±0.38, 
ax = 1.39±0.26; (c) the energy injection sample, with ao,i = 0.25±0.12, «o ,2 = 1.26±0.26, 
ax ,i = 0.30 ±0.27, and «x ,2 = 1.35 ±0.24; (d) the jet break sample, with a 0 ,i = 0.77 ±0.18, 
ao ,2 = 1-66 ± 0.16, «x,i = 0.95 ± 0.16 and «x ,2 = 1-70 ± 0.19. 
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Fig. 9.— The observed temporal break change Aa = « 2 —«i compared against the theoretical 
value predicted from the closure relations. The square (black) and circle (red) data points 
correspond to X-ray and optical bands, respectively. 
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Fig. 10.— The distributions of the temporal break change A a = «2 — oi\ for various sub¬ 
samples in the Gold sample, and the A a and a distributions of the Silver sample. Best 
Gaussian fits are marked with dashed curves: (a) the Aa distributions of the entire Gold 
sample, with Aao = 0.94 ±0.23 and A ax = 0.88 ±0.28; (b) the A a distributions of different 
sub-samples in the Gold sample, with Aao = 1-05 ± 0.17 and Aax = 1-10 ± 0.21 for the 
energy injection breaks, and Aao = 0.75±0.22 and Aax = 0.75±0.22 for the jet breaks; (c) 
and (d) the Aa distributions for the BPL and SPL GRBs in the Silver sample, respectively. 
The data are dispersed. 
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Spectral index 


Spectral index 


Fig. 11.— The distributions of the spectral indices (3 of various sub-samples in the Gold 
sample and their best Gaussian fits (dashed curves): (a) the distributions for the entire 
Gold sample, with /3o = 0.70 ± 0.15, and f3x = 0.98 ± 0.15; (b) the /3 distributions of 
the GRBs with optical and X-ray bands constrained in the same spectral regime, with 
Po = 0.77 ± 0.19, [3x = 0.89 ± 0.15; (c) the (3 distributions for the GRBs with optical and 
X-ray bands constrained in different spectral regimes, with (3q = 0.68±0.18, fix — 1.01±0.14 
and A/5 = Px~Po = 0.37±0.18, which is consistent with the theoretical value 0.5 as expected 
for a cooling break; (d) the f3 distributions for energy injection, with Po — 0-78 ± 0.12, and 
Px = 1.01 ±0.13; (e) the P distributions for the jet break sample, with P Q = 0.59 ±0.11, and 
Px = 0.97 ± 0.08; (f) the P distributions for the SPL decay sample, with Pq = 0.74 ± 0.24, 
and Px = 0.95 ± 0.19; (g) the P distributions for the GRBs with an ISM medium, with 
Po = 0.72 ± 0.21, and Px = 0.98 ± 0.10; (h) the P distributions for the GRBs with a wind 
medium, with P Q = 0.70 ± 0.10, and Px = 1-00 ± 0.20. 
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Fig. 12.— The distribution of the inferred electron spectral index p from various sub-samples 
in the Gold sample, and their best Gaussian fits (dashed lines): (a) the entire Gold sample 
with p = 2.33 ± 0.48, and the sub-samples with the optical and X-ray bands constrained in 
the same (p = 2.58 ± 0.39) and different (p = 2.17 ± 0.44) spectral regimes, respectively; (b) 
the sub-samples with different light curve behaviors: energy injection (p = 2.34 ± 0.38); jet 
break (p = 1.91 ±0.37); and SPL decay (p = 2.48 ±0.47); (c) the sub-samples with different 
medium types: ISM ( p = 2.43 ± 0.57) and wind ( p = 2.28 ± 0.33). 
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log (t b /ks) 


Fig. 13.— The distributions of the inferred % for the entire Gold sample, the energy injection 
sample, and the jet break sample. The dash lines are the best Gaussian fits: the entire Gold 
sample (log(tb/ks) = (3.8±0.9)), the energy injection break sample (log(tb/ks) = (3.6±1.9)) 
and the jet break sample (log(tb/ks) = (3.9 ± 0.7)). 
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q 


Fig. 14.— The distributions of the inferred energy injection parameter q and their best 
Gaussian fits (dashed lines): the entire Gold sample (q = 0.22 ± 0.11), the ISM sample 
(q = 0.20 ± 0.12) and the wind sample (q = 0.23 ± 0.13). 
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Fig. 15.— The derived upper limits of es of GRBs that are in spectral regime If {y m <v< 
u c ). Other parameters are fixed as e e = 0.1, n — 1, and A* = 1. The dashed line denotes 
6b = 10~ 5 . Most GRBs are consistent with this value, with a few having upper limits even 
below this value. 
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Fig. 16.— The distributions of various energy components derived from the Gold sample and 
their best Gaussian fits (dashed lines): (a) isotropic 7 -ray energy, E lr iso , with a typical value 
log(i7 7 !iso /erg) = (53.15 ± 0.69); (b) the total isotropic kinetic energy at the end of energy 
injection for the energy injection sample, with a typical value log (Ex, end/erg) = (54.99±0.86); 
(c) the distribution of the isotropic injected energy in the energy injection sample, with 
log(i?K,inj/erg) = (54.95 ± 0.61); (d) the isotropic kinetic energy at the deceleration time 
for the energy injection sample (log(£ , K,dec/erg) = (53.29 ± 0.45)), and for the entire Gold 
sample (log(77 K ,dec/ er g) = (54.66 ± 1.18)). The following parameters are adopted in the 
kinetic energy calculations: e e = 0.1, n = 1 or A* = 1, Y = 1, and eg = 10 -5 . 
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Fig. 17.— The correlations among different energy components: (a) E Kend — -E K ,inj : 

-EK,inj,52 = 0-69E , K i end°52 2 ) O 3 ) -^K,dec _ -ElCdiij: -EK,inj,52 = 41.777^^,52°) ( C ) ^ 1 ,iso — -EK,end : 

-EK,end,52 = 476.3; ?^52’ 23 ; (d), (e) and (f): E 7 , iso - E Kjdec for the entire Gold sam¬ 


ple, the energy injection sub-sample, and the sub-sample without energy injection, with 


Ei 


K,dec,52 


= 56.2 E\ 


-0.93±0.19 
7,iso,52 ? 


E] 


K,dec,52 


_ o q/?l-10±0.29 

— 0.»n/ 7iso52 


and E 


K, dec, 52 


= 316.2£™ 2 21 , respec¬ 


tively. 
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Fig. 18.— The radiative efficiency r\^ of the Gold sample as a function of E liiso (upper row) 
and histograms (lower row): (a) and (b): For the entire Gold sample. Among them, the 
GRBs in the energy-injection sub-sample are marked in red and black. Two efficiencies are 
calculated for each GRB: one using R K; dec (red) and the other using E end (black). Those 
GRBs without energy injection are marked in blue. Log-normal fits to the efficiencies derived 
from -Ek.dec and E end give log(?/ 7i d ec /%) = 0.75 ± 0.86 and log(r] 7i end/%) = —0.89 ± 0.97, 
respectively; (c) and (d): for the sub-sample of GRBs without energy injection, with best fit 
log(r/ 7 dec /%) = 0.36 ± 0.61. Open circles denote the ISM medium and solid squares denote 
the wind medium; (e) and (f): for the sub-sample of GRBs with energy injection, with best 
fits log(?7 7] dec/%) = 0.92 ± 1.25 and log(77 7i end/%) = —0.89 ± 0.97, respectively. Again ISM 
and wind cases are denoted as open circles and solid squares, respectively. 
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Fig. 19.— The distributions of the jet opening angle (9j), geometrically corrected 7 -ray 
energy (E y ) and kinetic energy (E K ), respectively, derived from the jet break sub-sample 
of the Gold sample. The dashed lines are the best Gaussian fits, with 6j = (2.8 ± 1.5)°, 
log(.E 7 /erg) = (49.86 ± 0.65), and log(E K /erg) = (50.89 ± 0.54), respectively. 























